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SECTION  I 
INTRODUCTION 


1.  1  OBJECTIVE 

The  objective  of  this  work  is  the  experimental  determination  of 
gas  properties  at  high  temperatures  and/ or  pressures.  The  tempera¬ 
ture  range  is  up  to  1500  K  and  the  pressure  range  is  up  to  400  atm. 

These  measurements  are  needed  for  engineering  needs  as  well  as  to 
support  potential  function  and  dense  gas  theory  requirements. 

Transport  and  thermodynamic  properties  are  the  data  base  for 
aerodynamics,  gas  phase  engine  design,  refrigeration,  and  many  im¬ 
portant  industrial  applications.  Uncertainties  of  15%  are  apparent!  in 
the  specific  heat  of  a  simple  gas  such  as  argon  at  ~800  atm  300  C.  The 
thermal  conductivity  of  many  common  polar  and  polyatomic  gases  and 
their  mixtures  are  not  known  to  50%  at  extremely  low  (100  C)  and  high 
temperatures  or  at  high  densities.  In  addition,  new  propulsion  systems 
and  laser  systems  generate  new  gas  property  requirements.  Recent 
experience  in  chemical  laser  research,  for  instance,  has  generated 
property  requirements  in  gases  for  which  no  information  was  available. 
Finally,  improvements  in  computational  and  test  capabilities  continually 
require  extended  range  or  accuracy  in  gas  property  information.  The 
long  lead  time  required  to  produce  new  measurements  and  to  critically 
evaluate  them,  calls  for  a  basic  research  effort. 

New  techniques  based  on  acoustic  and  schlieren  interferometric 
measurements  have  been  developed  at  Panametrics  to  help  expand  this 
data  base.  An  ultrasonic  pulse^»  3  technique  in  arcs  and  shock  tubes 
was  used  to  measure  transport  properties  of  largely  dissociated  and 
monatomic  gases.  The  ultrasonic  technique  becomes  intractable  in  the 
temperature  region  corresponding  to  about  0.  1%  to  80%  dissociation. 
Therefore,  a  schlieren  interferometric  technique^'  ®  employing  the  end 
wall  boundary  layer  in  a  shock  tube  was  developed. 

The  acoustic  pulse  technique  was  used  to  help  fulfill  the  require¬ 
ments  for  high  pressure  transport  and  thermodynamic  properties.  It 
was  found  that  the  sound  speeds  measured  this  way  could  provide  useful 
specific  heat  data  (~0.  1%).  However,  the  sound  absorption  consisted  of 
a  sum  of  viscosity,  thermal  conductivity  and  bulk  viscosity.  The  bulk 
viscosity  was  approximately  as  large  as  the  density  effect  in  viscosity 
and  thermal  conductivity,  therefore,  a  more  complete  theory  than  is 
available  today  is  required  to  make  pulse  absorption  measurement  useful. 
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However,  the  pulse  technique  can  provide  useful  information  on  dif¬ 
fusion  of  helium  and  hydrogen  in  the  heavier  gases.®’  ? 

The  data  required  to  satisfy  the  high  pressure  requirements 
must  be  high  precision  so  that  an  accurate  estimate  of  the  excess 
transport  properties  can  be  obtained.  Also,  independent  measure¬ 
ments  of  viscosity  and  thermal  conductivity  are  required.  These  ob¬ 
jectives  can  be  in  principle  achieved  with  an  acoustic  cavity  resonance 
technique.  High  precision  can  be  achieved  with  the  cavity  technique 
because  all  measurements  can  be  reduced  to  frequency.  The  center 
frequency  provides  sound  speed,  which  can  be  used  at  low  pressures 
to  calibrate  temperature  sensors.  The  linewidth  provides  transport 
properties.  Measurement  of  linewidth  for  different  types  of  reso¬ 
nances,  i.  e.  ,  radial,  longitudinal,  provides  simultaneous  viscosity 
and  thermal  conductivity  values.  Thus,  one  obtains  a  complete  set 
of  highly  accurate  thermodynamic  and  transport  properties  from  this 
cavity  technique. 

During  the  past  three  years  the  acoustic  resonant  cavity  tech¬ 
nique  for  measuring  the  thermodynamic  and  transport  properties  of 
gases  in  the  "low"  temperature  range  300  <  T  <  1000  K  and  pressure 
range  1-500  atm  has  been  demonstrated.  This  technique  compliments 
the  pulse  ultrasonic  optical  techniques  developed  for  elevated  tempera¬ 
tures  (1000  <  T  <  20,  000)  during  the  first  few  years  of  the  present 
effort.  The  cavity  technique  is  extremely  precise  because  the  mea¬ 
surements  have  been  reduced  to  easily  measured  quantities  -  frequency, 
a  90  phase  angle,  and  cavity  radius.  The  theory  of  the  resonator  in¬ 
dicates  that  the  acoustic  resonance  technique  should  be  capable  of 
yielding  0.  01%  in  viscosity  and  0.  05%  in  thermal  conductivity.  How¬ 
ever,  we  have  only  pushed  the  precision  to  0.  1%  in  viscosity  and  0.  5% 
in  thermal  conductivity  for  the  sake  of  supplying  immediate  engineering 
needs.  It  should  also  be  noted  that  0.  5%  viscosity  and  1%  thermal  con¬ 
ductivity  measurements  seem  possible®  with  wave  guide  driven  cavities 
up  to  ~2000°K. 

The  present  report  discusses  the  precision  of  the  cavity  tech¬ 
nique,  some  experimental  realizations  of  the  cavity  technique  electronic 
readout  options,  and  the  results  of  a  systematic  measurements  program 
in  nitrogen,  air,  hydrogen,  and  argon.  The  first  section  of  this  report 
discusses  how  well  the  cavity  response  can  be  described.  This  includes 
the  dependence  of  linewidth  and  center  frequencies  on  viscosity,  thermal 
conductivity,  relaxation  times  and  sound  speed  (thermodynamics).  Also 
discussed  here  are  the  effects  of  nonideal  boundary  conditions  such  as 
slip,  accommodation,  and  wave  leakage  into  the  cavity  walls.  The 
effects  of  the  transducers  on  the  cavity  response  are  discussed  in  the 
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second  section.  Mode  coupling,  transducer  coupling,  electrical  and 
acoustic  feedthrough,  are  included  in  the  cavity  response.  Tempera¬ 
ture  limitations  imposed  mainly  by  the  transducer  are  discussed.  All 
these  effects  are  dependent  on  how  one  chooses  to  construct  and  drive 
the  cavity.  Methods  of  "reading  out"  sufficient  information  electron¬ 
ically  to  determine  transport  and  thermodynamic  properties  are  given 
with  typical  results.  Finally,  the  experimental  gas  properties  deter¬ 
mined  in  the  past  year  and  a  half  are  presented. 
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SECTION  II 

THEORY  OF  IDEAL,  CAVITY  RESPONSE 

2.  1  LINEWIDTH  AND  CENTER  FREQUENCY  IN  ACOUSTIC 
GAS  RESONATORS 

The  linewidth  and  center  frequency  of  an  acoustic  resonance 
line  illustrated  in  Fig.  1  can  be  related  to  the  properties  of  the  gas 
in  the  cavity.  Since  we  are  using  small  amplitude  disturbances 
(1  pt  in  10^  to  10^)  the  linearized  conservation  equations  of  fluid 
dynamics  describe  the  cavity  performance  to  the  limit  that  the  con¬ 
tinuum  theory  is  valid.  In  addition,  the  linearity  allows  the  device 
equations  to  be  worked  out  as  accurately  as  one  desires.  First,  the 
completeness  of  the  device  equations  is  never  in  question  as  is  the 
case  with  the  nonlinear  corrections  (convection,  edge,  and  entrance 
corrections)  necessary  in  other  techniques  for  transport  property 
measurement.  Second,  the  calculations  can  always  be  carried  out 
with  orthogonal  functions  and  perturbation  theory  in  such  a  way  that 
the  maximum  error  in  the  neglected  terms  can  be  calculated.  Finally, 
slip  and  accommodation  can  easily  be  included  in  the  theory  by  mod¬ 
ifying  the  wall  boundary  conditions.  The  section  below  will  show  how 
the  linewidth  and  center  frequency  are  related  to  the  gas  properties 
and  boundary  conditions  to  0.  01%  in  linewidth.*  However,  before 
launching  into  details,  a  summary  of  the  first  order  results  is 
appropriate. 

The  gas  property  information  is  contained  in  the  resonator  line 
shape.  The  line  shape  is  the  complex  pressure  amplitude  versus  fre¬ 
quency  characteristics  of  the  cavity.  Although  the  line  shape  can  have 
a  very  complex  mathematical  form,  a  Lorentzian  line  shape  is  a  useful 
approximation  because  of  its  simplicity.  The  pressure p  in  the  sound 
field  of  a  cavity  driven  at  frequency  f  is 

-  _  _ AC  . 

P'2ir(f-f)  +  i5 
o 

in  the  Lorentzian  approximation.  A  is  the  pressure  at  the  maximum 
amplitude  point,  that  is,  at  f  =  fQ.  The  center  frequency  is  fQ,  and 
the  half  width  is  6  . 


*This  corresponds  to 


0.  01% 
Q 


in  sound  speed. 


See  below. 
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Fig.  1.  Variation  of  phase  and  amplitude  with  frequency  near  a  resonance. 
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The  complex  pressure  is  used  to  simplify  the  bookkeeping,  and 
is  composed  of  the  pressure  amplitude  and  a  phase  factor  ex9  .  Only 
complex  pressures  are  found  in  nature  because  phase  in  addition  to 
frequency  and  pressure  is  necessary  to  specify  the  sound  wave.  The 
phase  of  the  complex  pressure  in  the  cavity  is  measured  relative  to 
the  pressure  wave  driving  the  cavity.  Experimentally,  the  cavity  is 
driven  by  some  incident  pressure  wave  at  one  end  wall  and  the  result¬ 
ing  field  in  the  cavity  sensed  at  the  other  end  wall.  The  frequency  of 
the  driving  wave  is  varied  while  the  amplitude  and  phase  of  the  pres¬ 
sure  in  the  cavity  are  measured. 


The  amplitude  and  phase  variation  versus  frequency  is  shown 
in  Fig.  1  near  an  acoustic  resonance.  The  most  distinctive  feature 
of  the  cavity  response  is  the  maximum  amplitude  at  fQ.  The  center 
frequency  can  be  measured  using  this  amplitude  measurement  to  a 

j  which  depends  on  the  linewidth  (5  )  and  the  amplitude  uncertainty 


by  the  relation 


Af 
_ _ o 

f 

o 


6  dA 

2  TT  f  A 
o 


(2) 


where  A  A  is  the  pressure  amplitude  jitter  and  A  is  the  pressure 
amplitude.  The  linewidth  can  be  measured  by  varying  the  frequency 
on  either  side  of  the  center  frequency  until  the  amplitude  reaches 

— ■=:  A  as  indicated  in  Fig.  1 . 

Similar  definitions  are  possible  in  terms  of  the  phase  versus 
frequency  curves.  The  half  power  points  correspond,  for  a  Lorentzian 
line,  to  the  45°  and  135°  phase  shift  points.  In  fact,  the  phase  is  ulti¬ 
mately  used  in  the  present  work  to  define  the  linewidth  because  the  phase 
measurement  can  be  made  more  accurately  than  the  amplitude  determi¬ 
nation.  This  is  discussed  below  in  the  section  on  experimental  technique. 

The  cavity  has  several  types  of  resonances  or  modes,  each  being 
an  independent  solution  to  the  equations  of  motion  and  the  cavity  boundary 
conditions.  The  relative  contributions  of  viscosity  and  thermal  conduc¬ 
tivity  to  the  linewidth  are  different  for  each  mode.  Thus,  measurement 
of  two  modes  allows  the  transport  properties  to  be  separately  determined. 
The  cavity  dimensions  determine  the  frequency  at  which  each  mode  is 
appreciably  excited.  Thus,  as  we  increase  frequency  the  various  modes 
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rise  to  a  maximum  and  fall  to  zero  within  a  few  linewidths  of  their 
center  frequency.  A  typical  amplitude  versus  frequency  plot  is  shown 
in  Fig.  2.  Any  mode  can  be  isolated  by  driving  the  cavity  near  the 
resonant  frequency  of  the  mode  in  question. 


The  first  order  expressions  for  sound  speed  and  linewidth  in 

terms  of  the  gas  properties  are  accurate  to  the  order  of 
The  center  frequency  is 


f.. 

k) 


a 


c 

2d 


The  linewidth  is 


6.. 


f.. 

ij 


1/2 


6  f.. 
v  ij 


1/2 


(3) 


(4) 


Here  the  first  subscript  refers  to  the  radial  mode  number  and  j  refers 
to  the  longitudinal  mode  number.*  The  c  is  the  sound  speed,  d  the  length 

of  the  cavity,  j  the  kinematic  viscosity,  y  the  specific  heat  ratio, 

Pr  the  Prandtl  number  and  6^.  the  volume  loss.  The  quantitatives,  a-, 
ajj  and  b^j,  are  known  quantities  which  depend  on  cavity  geometry  and  J 
the  mode  numbers.  These  constants  as  well  as  5y  are  given  explicitly 
in  Table  1.  Typical  linewidth  and  center  frequency  values  are  120  Hz  in 
20  kHz  at  1  atm,  decreasing  to  1  5  Hz  in  30  kHz  at  400  atm. 

Equations  (3)  and  (4)  describe  the  losses  due  to  dissipation  in  the 
cavity  itself.  Additional  losses  discussed  in  later  sections  occur  due  to 
the  transducers  and  mode  coupling.  The  cavity  is  designed  so  that  the 
dissipation  in  the  cavity  itself  is  the  dominant  loss.  The  rest  of  this 
section  is  devoted  to  the  accuracy  of  these  two  relationships. 


*(jo)  are  the  radial  modes,  (oj)  are  the  longitudinal  modes,  and  (ij)  are 
compound  modes.  See  Jack  Bradshaw's  thesis  9  and  Ref.  3  for  more 
details. 
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Table  1.  Constants  for  first  approximation  linewidths  for  symmetric 
modes  in  a  cylindrical  acoustic  cavity. 
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2.  2  PHYSICAL  BASIS  OF  CAVITY  EQUATIONS 

Following  Thiesen*^  and  Fritsche** '  the  exact  equivalence 
of  the  conservation  equations  of  fluid  dynamics  with  a  set  of  three  wave 
equations  was  demonstrated  in  a  previous  report. ^  These  three  wave 
equations  describe  the  propagation  of  acoustic,  viscous,  and  thermal 

lie  “ 

waves.  These  waves  propagate  independently  of  each  other  except  at 
finite  amplitudes  and  at  discontinuities  in  elastic  properties  (such  as  the 
cavity  walls).  The  boundary  conditions  on  the  cavity  walls,  continuity 
of  temperature  and  strain,  can  only  be  satisfied  if  thermal  and  viscous 
waves  are  generated  on  the  walls. 

The  boundary  conditions  on  the  walls  of  the  cavity  must  be 
carefully  evaluated.  Application  of  the  boundary  conditions  to  the 
three  wave  equations  is  then  only  a  formal  problem.  The  idealized 
boundary  conditions  used  by  Fritsche  assume  that  the  wall  has  much 
greater  capacity  to  absorb  heat  than  the  gas.  Thus,  the  temperature 
fluctuations  at  the  wall  vanish.  In  addition,  the  wall  is  assumed  to  be 
much  stiffer  than  the  gas  so  that  the  particle  velocity  vanishes  there. 
Due  to  velocity  slip,  thermal  accommodation,  finite  wall  thermal 
capacity  and  finite  wall  acoustic  impedance,  the  particle  velocity  and 
fluctuating  temperature  do  not  vanish.  These15,  ^  effects  have  been 
discussed  for  the  lowest  order  traveling  wave  in  a  tube.  Generaliza¬ 
tion  of  these  wall  effects  to  higher  order  traveling  waves  and  to  waves 
reflected  from  end  walls  is  given  in  Appendix  A  below.  However,  mean 
freepath  effects,  slip  and  accommodation  require  further  discussion 
here. 


The  wavelengths  X.  v  and  X  t  for  the  viscous  and  thermal  waves, 
respectively,  are  given  by3 


X 

v 


(5) 


The  continuum  theory  might  be  expected  to  be  in  error  in  an  amount 
roughly  equal  to  the  ratio  of  the  wavelength  to  the  mean  free  path 
squared  since  the  theory  includes  first  order  deviations  from  equilib¬ 
rium.  A  mean  free  path  (L)  can  be  defined  by  the  zero  pressure  vis¬ 
cosity  ti  and  the  expression 


*The  viscous  and  thermal  waves  at  the  vortex  and  entropy  are  used  in 
the  theory  of  turbulence.1  14  The  entropy  wave  is  also  responsible 
for  the  central  peak  of  Brilloiun  scattered  light. 
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L 


(6) 


Values  of  —  for  representative  conditions  (shown  in  Table  2)  indicate 

that  the  continuum  formulas  for  absorption  (=  linewidth)  should  be  good 
to  pt  in  1  0^. 

Slip^'  ***  and  accommodation^’  ^  are  mean  free  path  effects 
due  to  imperfect  momentum  and  energy  exchange  within  a  few  mean  free 
paths  of  the  wall.  If  the  molecules  were  simply  reflected  like  billiard 
balls  the  surface  would  look  like  a  stationary  layer  of  gas.  The  tempera* 
ture  of  the  wall,  for  instance,  would  simply  be  the  value  extrapolated 
from  the  temperature  gradient  in  the  wall.  However,  a  certain  fraction 
of  the  molecules  stick  to  the  wall  before  being  reemitted  into  the  hotter 
flowing  gas.  Thus,  the  gas  in  a  layer  a  few  mean  free  paths  thick  has  a 
lower  flow  velocity  and  a  higher  temperature’1'  than  one  would  predict 
from  the  gradients  by  extrapolation  to  the  wall.  The  size  of  this  effect 
depends  on  the  fraction  of  gas  molecules  which  become  attached  to  the 
wall.  This  is  mathematically  equivalent  to  (but  not  physically  equal  to) 
assigning  a  length  i  y  and  f  t  for  the  reflected  molecules  to  adjust  to  the 
flow  and  temperature  of  the  gas. 


The  effect  of  slip  is  to  give  a  net  flow  velocity  on  the  wall, 
wall  flow  with  slip  is 


The 


U 


3'U 
v  9  r 


(7) 


instead  of  U  =  0.  Similarly,  the  fluctuating  component  of  temperature 
at  the  wall  is 


T 


9  T 
a  r 


(8) 


instead  of  zero.  The  correction  to  the  linewidth  which  results  from  slip 
and  accommodation  is  the  same  as  that  in  Table  1  with  l  and  l  ^  in  place 
of  mean  free  path.  Since  1  v  and  i  ^  are  a  few  mean  free  paths  typically 
(and  can  be  as  much  as  10  L)  these  corrections  are  somewhat  greater  than 
those  given  in  Table  2.  The  measurements  reported  below  were  taken 


*The  temperature  is  somewhat  lower  for  a  negative  temperature  gradient. 
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Table  2.  Ratio  of  mean  free  path  to  viscous  wavelength  for 


nitrogen  under  typical  pressure-temperature  and 
frequency  conditions. 


j  PRESSURE  ATM  | 

0.  1 

10  kc 

Temperature 

1  00  kc 

Temperature 

300 

1000 

300 

1000 

2.  3  x  10  2 

4.  92  x  10_1 

4.  00  x  10  2 

8.  52  x  10”1 

10 

2.  20  x  10"3 

2/11  x 10"3 

3.  81  x  10"3 

3.  65  x  10  3 

500 

6.  47  x  1  0  k 

1.  71  x  10“5 

1.12x10  3 

2.  96  x  1 0 
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under  frequency,  pressure  and  temperature  conditions  where  slip  and 
accommodation  coefficients  contribute  less  than  1  pt  to  10^  to  the  vis¬ 
cosity  and  thermal  conductivity.  However,  under  the  as  sumption  that 
the  theory  for  slip  and  accommodation  corrections  given  in  Appendix  A 
is  good  to  10%,  it  is  possible  to  correct  measurements  of  transport 
properties  to  1  pt  in  10^  for  densities  between  10-100  amagats. 


The  losses  to  the  side  walls  and  end  walls  due  to  acoustic  trans¬ 
mission  into  the  cavity  walls  give  rise  to  a  linewidth  component  Anp 
which  must  be  added  to  Eq.  {4)  to  get  the  total  cavity  linewidth.  It  is 
difficult  to  calculate  this  effect  exactly  since  the  type  of  acoustic  motion 
excited  in  the  cavity  walls  can  be  very  complex.  Nevertheless,  as  shown 
in  Appendix  A,  the  linewidth  due  to  acoustic  transmission  is  proportional 
to  frequency  and  the  ratio  of  the  acoustic  impedance  of  the  gas  (p  c)  to  the 
acoustic  impedance  of  the  cavity  walls  (p^c^). 

The  acoustic  transmission  linewidth  is 


A 

np 


a  f  + 
s 


a . 


P  c 
w  w 


(9) 


where  a  g  and  «£  are  constants  for  each  cavity.'1'  In  addition,  a  s  and  a ^ 

are  approximately  equal  to  for  the  present  cavities.  The  acoustic 

impedance  ratio  varies  from  ~10"^  at  1  atm  to  10  ^  at  100  atm  for  the 
steel  cavities  used  in  the  present  work.  Note  the  acoustic  transmission 
losses  to  the  end  wall  goes  to  zero  for  a  radial  resonance.  This  is  be¬ 
cause  the  wave  number  kpU)  goes  to  zero  at  a  radial  resonance  and  there 
is  no  axial  particle  motion  to  couple  to  the  wall. 

In  practice,  a  g  is  determined  by  measuring  more  than  one  radial 
mode  in  the  cavity  of  interest.  The  difference  between  the  total  linewidth/ 
(frequency)l then  provides  a  s.  Similarly,  a  £  can  be  measured  by  the 
same  technique  using  more  than  one  longitudinal  mode.  Notice  that  a  g 
and  <*■£  do  not  depend  on  the  gas  so  that  a  self-consistent  check  is  obtained 
from  measurements  in  different  gases. 

The  linewidth  component  Anp  can  be  obtained  to  ~0.  5%  in  the  present 
experiment.  Thus,  an  uncertainty  aue  to  Anp  of  typically  1  pt  in  10^  in  vis¬ 
cosity  at  1  atm  becomes  1  pt  in  10^  at  100  atm  and  increases  to  3  pts  in  10^ 

*This  assumes  that  there  are  no  resonances  in  the  cavity  walls  at  the 
operating  frequencies. 
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at  500  atm.  There  are  several  ways  to  reduce  the  size  of  the  acoustic 
transmission  loss  which  we  have  not  had  time  to  realize.  For  instance, 
high  density  and  modulus  material  such  as  tungsten  can  be  used  for  cavity 
construction.  In  addition,  radial  compound  modes  and  especially  asym¬ 
metric  modes?  can  be  employed  instead  of  longitudinal  modes.  Finally, 
large  cavities  can  be  used  to  accentuate  the  transport  property  losses 


(which  vary  as 
1  . 

vary  as  —  ). 


^==-  )  compared  to  the  acoustic  transmission  losses  (which 


The  losses  calculated  from  the  linewidth  expressions  given  above 
are  somewhat  larger  than  the  real  dissipation  in  the  cavity.  The  wave 
equation  cannot  be  separated  exactly  into  radial  wave  functions  due  to  the 
corners  of  the  cavity.  For  instance,  when  a  traveling  wave  which  satis¬ 
fies  the  side  wall  boundary  conditions  exactly  is  incident  on  the  end  wall 
a  new  set  of  viscous  and  thermal  waves  is  required  (A.  56)  to  satisfy  the 
end  wall  conditions.  These  waves  are  assumed  to  propagate  10“  ^  cm  or 
less  away  from  the  cavity  wall.  The  waves  originating  on  the  end  wall 
lying  near  the  corner  give  rise  to  unbalanced  temperature  and  normal 
velocity  on  the  side  wall.  Similarly,  the  viscous  and  thermal  waves 
originating  from  those  portions  of  the  side  wall  near  an  end  wall  (A.  50) 
cause  similar  imbalanced  temperatures  and  velocities  on  the  end  wall. 


Expansions  of  cavity  wave  functions  can  be  used  to  satisfy  the 
boundary  conditions  more  exactly.  Radial  wave  functions  may  be  used 
to  satisfy  the  end  wall  boundary  conditions  and  longitudinal  mode  wave 
functions  may  be  used  to  provide  a  suitable  expansion  from  the  side  wall 
boundary  conditions.  A  similar  problem  arises  for  acoustic  resonances 
in  a  finite  cylindrical  solid  resonator  when  one  tries  to  satisfy  shear 
and  normal  stress  boundary  conditions. »  ^2 


A  good  approximation  can  be  obtained  by  noting  that  the  corners 
must  be  temperature  and  velocity  nodes.  In  addition,  the  troublesome 

-k.r  -k.z 

waves  vanish  roughly  as  e  or  e  ,  (i  =  2,  3),  for  the  end  wall 

and  side  wall  misfits,  respectively.  Thus,  we  can  assume  that  in  the 
final  solution  the  energy  density  (E)  must  be  a  function  of  distance  along 
the  wall  instead  of  constant  (as  assumed  in  Appendix  A).  On  the  side 

-|  k^z | 

wall  the  viscous  wave  amplitude  must  fall  off  as  e  with  distance 

from  the  end  wall.  Similar  relations  hold  for  the  thermal  waves  and  on 
the  end  walls  for  both  types  of  waves.  Integrating  the  reduced  energy 
density  over  the  wall  area  a  correction  factor  for  viscosity 
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and  thermal  conductivity 


are  obtained.  The  correction  factors  were  always  less  than  1  pt  in  10 
in  the  present  work  and  were  neglected. 
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SECTION  in 

CORRECTIONS  TO  BASIC  CAVITY  THEORY  DUE  TO  TRANSDUCER 
COUPLING  AND  ELECTRICAL  FEEDTHROUGH 

3.  1  EFFECTS  OF  COUPLING  HOLE  TRANSMISSION  LINES 

The  linewidth  of  a  resonance  is  a  measure  of  how  fast  energy 
leaks  out  of  the  standing  wave  trapped  in  the  acoustic  cavity.  The 
total  linewidth  consists  of  the  "natural"  linewidth  due  to  dissipation 
of  energy  through  thermal  conduction  and  viscous  friction  (discussed 
in  Section  II  plus  additional  linewidth  due  to  energy  transmitted  out 
of  the  cavity.  For  instance,  in  a  completely  closed  cavity  the  trans¬ 
mission  loss  to  the  finite  acoustic  impedance  of  the  end  wall  is  ~0.  01 
cycles  at  1  atm  and  increases  to  1  cycle  at  100  atm. 

Similarly,  when  the  transducers  are  used  to  couple  sound  in 
and  out  of  the  cavity  there  are  additional  linewidth  components  due 
to  energy  which  propagate  into  the  transducers.  The  phase  shift  due 
to  the  transit  of  sound  from  the  electrical  drive  to  the  cavity  on  the 
transmit  side  and  a  similar  phase  shift  on  the  receive  side  add  to  the 
resonance  phase  shift  in  the  signal.  In  addition,  the  electrical  signal 
used  to  drive  the  transducers  can  take  a  path  through  the  cavity  walls 
or  directly  by  electromagnetic  pickup.  We  measure  the  phase  versus 
frequency  to  obtain  linewidth.  The  phase  shift  due  to  each  of  these 
effects  give  rise  to  an  effective  linewidth  component  which  must  be 
processed  out  to  obtain  transport  properties. 

The  relative  size  and  sign  of  each  of  the  above  linewidth  com¬ 
ponents  are  displayed  in  Table  3  for  the  type  of  cavities  employed  in 
the  present  experiments.  An  examination  of  Table  3  indicates  that  the 
hole  size  of  0.  004"  or  smaller  is  desirable'*  at  elevated  pressures.  In 
addition,  accurate  calculation  of  correction  factors  due  to  the  hole  is 
necessary  to  provide  accurate  absolute  viscosity  and  thermal  conductivity 
measurements.  The  sound  was  coupled  into  the  cavity  through  small 
holes  in  the  present  work.  Although  the  corrections  for  the  hole  in  the 
cavity  are  difficult  to  calculate  (for  a  0.  1-0.  01%  linewidth  measurement), 
this  coupling  technique  has  many  practical  advantages  over  other  possible 
transducers. 


^However,  hole  sizes  greater  than  0.  008"  were  used  in  the  present  work. 
Future  work  should  include  the  development  of  smaller  coupling  hole 
transducers. 
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Table  3.  Typical  linewidth  components  (argon  300°K)  for 
the  first  longitudinal  mode. 


Pressure 

Hole 

diameter 

Natural 

linewidth 

Hole  loss 
linewidth 

Transducer 
phase  shift 

6 

5 

6 . 

o 

a 

h 

1 

0.  013 

119 

2 

0.  004 

0.  05 

0.  5 

100 

0.  013 

12.  5 

27 

2 

0.  004 

0.  1 

0.  3 

500 

0.  013 

9.5 

27 

2 

0.  004 

2 

0.  3 
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The  losses  to  the  transducer  in  the  present  case  consist  of 
several  terms.  In  general,  the  acoustic  impedance  of  the  coupling 
hole  terminated  by  the  electrostatic  transducer  could  be  calculated 
in  the  usual  way.^3,  24  T^g  reflection  coefficient  so  generated 
could  be  put  into  the  usual  formulas  for  linewidth  of  a  transmission 
line  (which  the  cavity  is)  terminated  by  the  impedance  of  coupling 
holes.  However,  we  have  found  detailed  tracing  of  multiple  reflec¬ 
tions  much  simpler  to  use  in  the  present  case.  This  is  particularly 
true  when  it  comes  time  to  see  which  physical  parameters  are  sig¬ 
nificant  to  sort  out  approximations  and  to  calculate  second  order 
terms. 


The  increase  in  linewidth  due  to  the  traveling  wave  multiple 
reflections  depends  only  on  the  complex  reflection  coefficients  (R) 
calculated  in  Appendix  C.  These  reflection  coefficients  include- 
the  loss  of  energy  due  to  the  excitation  of  other  modes  during  each 
reflection.  The  contribution  of  these  other  modes  is  considered  in 
Section  3.  3.  The  coupled  modes  behave  like  feedthrough  provided 
the  cavity  modes  are  separated  by  10  linewidths.  Finally,  the  addi¬ 
tional  phase  shift  due  to  multiple  reflections  in  the  coupling  holes  is 
given  in  Appendix  B. 

The  total  signal  received  can  be  calculated  by  adding  up  the 
multiple  reflections  of  the  various  waves  trapped  in  each  of  the  three 
sections  of  the  cavity.  As  shown  in  Appendix  B  the  total  receive  sig¬ 
nal  is 

Cf  -  AoF!T1.  FT.2F2-  (10> 

Here  A0  is  the  amplitude  generated  by  the  vibration  of  the  drive  trans¬ 
ducer.  Fj  and  F 2  are  the  response  coefficients  which  result  from 
adding  up  all  multiple  reflections  in  the  transmit  hole  (F^ )  and  the 
receive  coupling  hole  (F2)  given  by  Eq.  (B.  1).  T^  and  T  2  are  the 
traveling  wave  transmission  coefficients  from  the  transmit  coupling 
hole  into  the  cavity  Tj  and  from  the  cavity  into  the  coupling  hole  on 
the  transmit  side  Tt  2*  The  reflection  coefficients,  Rj^  and  R22»  as 
well  as  Tj  and  T  2  (see  Eq.  (B.  1)),  are  calculated  from  the  equations 
in  Appendix  C.  The  additional  reflection  coefficients,  Rqq  and  R33, 
acquired  to  calculate  F^  and  F2#  respectively,  are  taken  as  unity.  More 
exact  values  for  these  reflection  coefficients  have  been  calculated  from 
the  theory  of  electrostatic  transducers. ^5  Measurements  of  linewidth 
with  the  transducers  reversed  were  made,  since  an  undefined  constant 
is  present  in  the  transducer  theory. 
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Only  the  phase  of  the  quantity  T  £  F  2  is  needed  to 

calculate  the  effect  of  phase  shift  in  the  coupling  hole.  This  phase 
shift  is  always  less  than  ~25  in  the  present  experiment  because  the 
first  resonance  in  Fj  or  F2  occurs  at  3  to  10  times  the  operating 
frequency  for  the  cavity  modes  employed  in  the  present  measurement. 
The  change  in  coupling  hole  phase  shift  (0^)  is  always  less  than  2.  5 
over  a  typical  linewidth.  This  is  because  the  Q  of  the  coupling  hole  is 
always  less  than  one -tenth  of  the  cavity  Q.  The  effective  linewidth  of 
this  correction  is 


6 


h 


(ID 


This  is  the  quantity  listed  under  transducer  phase  shift  correction  in 
Table  3. 


3.  2  EFFECTS  OF  COUPLING  HOLE  REFLECTION 
COEFFICIENTS 


The  largest  correction  due  to  the  coupling  holes  is  contained  in 
(F)  itself.  Considering  the  end  walls  perfectly  reflecting  except  for 
the  losses  to  the  coupling  holes  and  mode  coupling,  the  linewidth  due 
to  the  coupling  hole  can  be  written 
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2  f 
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nir 


1  -  r 

_ EE 

2 

1  +  r  pn 


(12) 


where  fQn  is  the  resonant  frequency  of  the  nth  longitudinal  mode  and  r 
is  the  reflection  coefficient  amplitude  from  Appendix  C.  This  linewidth 
may  be  simply  added  to  the  other  linewidth  components  described  above. 


The  linewidth  due  to  the  coupling  hole  represented  by 


6 

a 


w 


(13) 


where  G  (pn)  is  a  constant  for  each  mode,  rQ  is  the  coupling  hole  radius, 
and  a  the  cavity  radius.  The  loss  component  6a  represents  the  energy 
simply  propagated  into  the  hole.^6 

A  second  component  of  the  5pn  arises  from  the  complex  parts  of 
the  wave  numbers  as  described  in  Appendix  C.  This  linewidth  component 
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gr7(P“)  +  gx(P“) 


(14) 


is  dependent  on  the  Prandtl  number  and  the  viscosity.  Here,  (pn) 
and  gx  (pn)  are  constants  for  each  mode  which  are  considered  to  be 
additional  factors  in  the  desired  cavity  response.  Thus,  the  hole 
correction  is  only  5a  (see  below). 

3.  3  MODE  COUPLING  EFFECTS 


The  reflection  coefficients  of  Appendix  C  account  only  partially 
for  the  energy  lost  to  mode  coupling.  In  fact,  pressure,  carried  by 
these  other  modes  does  reach  the  receive  transducer.  However,  if 
the  modes  are  well  separated  the  additional  signal  due  to  mode  cou¬ 
pling  varies  only  slowly  and  looks  like  a  form  of  feedthrough. 

There  are  experimental  precautions  one  takes  to  avoid  mode 
coupling.  First,  the  overall  coupling  between  transducer  and  cavity 
is  reduced  as  much  as  possible  consistent  with  signal -to -noise  con¬ 
siderations.  Second,  the  symmetry  of  the  transducer  is  chosen  to 
favor  the  desired  modes.  For  instance,  an  electrostatic  transducer 
covering  the  entire  face  of  the  cavity  employed  in  the  present  work 
excites  the  first  20  longitudinal  modes,  preferentially  with  radial  and 
asymmetric  modes  down  by  at  least  40  dB  (a  factor  of  100).  Similarly, 
a  small  centrally  located  transducer  excites  longitudinal,  radial  and 
compound  modes,  the  asymmetric  mode  being  undetectable  (40  dB  down). 
Finally,  the  ratio  of  cavity  length  to  diameter  is  chosen  for  maximum 
frequency  separation  between  modes.  However,  after  all  precautions 
have  been  taken  one  is  typically  left  with  a  mode  coupling  phase  shift 
equivalent  to  ~0.  5  to  0.  05%  of  the  linewidth.  Thus,  to  realize  the  ac¬ 
curacy  inherent  in  the  cavity  technique  a  mode  coupling  correction  is 
necessary. 

The  response  of  a  cavity  can  be  calculated  in  great  detail  by 
considering  multiple  reflections  of  the  cylindrical  wave  guide  traveling 
waves  of  Appendix  A  from  the  cavity  end  walls.  Initially,  the  transducer 
excites  many  modes  according  to  the  transmitter  coupling  coefficient 
t^11  for  each  radial  mode  <J)n  as  indicated  in  Appendix  C.  When  the  fre¬ 
quency  of  the  driving  transducer  excitation  is  near  the  resonant  frequency 
of  the  sth  mode,  the  amplitude  of  the  resonant  mode  is  ~Qtms.  All  the 
other  modes  are  ~tmn  in  amplitude.  Each  time  the  wave  system  reflects 
from  one  of  the  end  walls,  reflection,  transmission  and  mode  conversion 
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occurs.  The  amount  of  each  wave  generated  depends  on  the  reflection 
coefficient  Rgs  and  the  coupling  coefficient  Rns  of  the  transducers 
involved. 


The  calculation  of  the  mode  coupled  signal  is  similar  to  the 
derivation  in  Appendix  B.  A  right  going  wave 


♦  .  = 


t 

m  m 


i(ut  +  k  ^z) 


(15) 


is  introduced  by  the  transmitter.  Here  0  is  the  radial  wave  function 
for  the  resonant  wave  in  question.  The  transmitted  wave  bounces  back 
and  forth  between  the  end  walls  of  the  cavity.  These  multiple  reflec¬ 
tions  give  rise  to  a  total  right  going  wave  for  the  resonant  mode  on  the 
receive  wall 
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t  tif 
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All  the  singly  mode  converted  waves  0n+  arising  from  the  reflection  of 
0  s  from  the  receive  transducer 


s 


waves  which  undergo  multiple  mode  conversions  are  proportional  to 
higher  order  powers  of  the  coupling  coefficients.  Since  the  coupling 
coefficients  are  small,  these  higher  order  terms  are  neglected.  A 
similar  contribution  of  singularly  mode  converted  waves  is  generated 
by  reflections  from  the  transmit  coupling  hole.  Thus,  the  total  mode 
converted  wave  source  strength  in  the  nth  radial  mode  is 

d >  =^tn+2^tSRnf  l  +  e 

r n  mm  m  m  s 

The  first  term  is  simple  transmission  into  the  cavity  from  the  transmit 
coupling  hole.  This  term  gives  rise  to  an  acoustic  feedthrough  which 

2 

contributes  less  than  of  the  total  received  signal  at  the  half  power 

points,  where  Q  is  the  cavity  quality  factor  j  •  This  feedthrough¬ 
like  component  is  processed  out  of  the  data  as  described  in  Section  3.  4. 
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The  second  term  in  Eq.  ( 1 6 )  is  not  mentioned  in  the  literature. 
However,  since  it  is  proportional  to  the  resonant  mode  response  func¬ 
tion  Fg,  it  can  dominate  the  mode  coupling  effects.  The  total  signal 
received  due  in  the  presence  of  mode  coupling*  is 
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The  second  term  in  brackets  represents  the  mode  coupling  contribution 
and  is  about  10"^  of  the  first  term.  When  the  modes  are  well  separated 
the  mode  coupling  term  does  not  depend  on  frequency  because  we  are  in 
the  wings  of  Fn.  Thus,  the  major  effect  of  the  mode  coupling  term  is  to 
add  a  system  phase  shift.  The  system  phase  shifts  are  removed  by  the 
iterative  determination  of  center  frequency  and  phase  shift  at  center 
frequency.  Nevertheless,  more  accurate  calculation  of  the  phase  shift 
due  to  mode  coupling  is  possible  using  the  coupling  and  transmission 
coefficients  from  Appendix  C,  although  this  was  not  done  with  the  data 
reported  below. 


3.  4  FEEDTHROUGH  ROUTINES  FOR  DATA  REDUCTION 


Summarizing  the  last  three  sections,  accurate  determination  of 
natural  line  shape,  that  is,  phase  and  amplitude  versus  frequency  for 
a  given  resonance,  requires  elimination  of  competing  phase  shift  and 
amplitude  variations.  The  effects  of  various  components  on  the  mea¬ 
sured  phase  and  amplitude  can  be  more  clearly  understood  by  reference 
to  a  diagram  of  the  experimental  apparatus  (Fig.  3).  A  frequency  syn¬ 
thesizer  is  used  to  generate  a  reference  signal  for  the  lock-in  amplifier 
and  a  driving  signal  for  the  resonant  cavity.  Between  the  synthesizer 
and  the  cavity,  a  line  amplifier  is  adjusted  to  give  enough  output  to  drive 
the  transmitted  transducer.  The  output  of  the  transmit  transducer  enters 
the  cavity  through  a  coupling  hole  and  stimulates  sound  waves.  The  am¬ 
plitude  of  the  sound  field  in  the  cavity  is  sensed  through  a  second  coupling 
hole,  and  is  picked  up  by  the  receive  transducer.  The  output  of  the  re¬ 
ceiver  is  fed  into  a  tuned  amplifier  to  increase  the  signal -to-noise  ratio. 
Finally,  the  lock-in  amplifier  provides  the  signals  necessary  to  determine 
the  amplitude  of  the  signal  and  the  phase  of  the  signal  with  respect  to  the 


'^Tj  has  been  used  in  place  of  tms  for  simplicity.  It  should  be  borne 
in  mind  that  Tj.  and  T^  2  depend  on  radial  mode  number.  For  the  modes 
used  in  the  present  work  all  modes  above  m  =  0  in  the  coupling  hole  are 
below  cutoff  and  contribute  a  negligible  amount  of  signal. 


36 


OJ 


Fig.  3.  Signal  paths  for  resonant  cavity  transport  and  thermodynamic  property  measurements. 
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electrical  driving  voltage.  The  amplitude  information  is  read  out  on 
a  digital  voltmeter.  The  phase  is  measured  either  on  a  phase  shift 
dial  in  the  reference  channel  or  preferably  as  in-phase  and  out  of 
phase  components  of  amplitude  on  a  digital  voltmeter. 

The  major  objective  of  the  data  reduction  discussed  here  is  the 
elimination  of  the  systematic  phase  shift  and  feedthrough.  The  phase 
shift  due  to  all  signal  components  other  than  those  due  to  the  resonant 
cavity  has  to  be  eliminated  in  order  to  determine  the  true  cavity  line- 
width  and  line  shape.  Systematic  phase  shift  is  introduced  by  the 
following  elements  of  the  apparatus:  The  frequency  synthesizer,  line 
amplifier,  coupling  holes,  tuned  amplifier,  and  distributed  RC1  s.  All 
these  component  phase  shifts  are  in  series  and  add  up  numerically. 

In  addition,  the  amplitude  variations  with  frequency  due  to  elements 
in  the  signal  path  can  be  corrected  by  multiplication  by  a  frequency 
dependent  numerical  factor. 

A  second  set  of  interfering  signals  are  the  acoustic  and  elec¬ 
trical  feedthroughs  "around"  the  cavity.  The  feedthrough  signal  path 
is  parallel  to  the  signal  path  indicated  in  Fig.  3  by  the  imaginary  re¬ 
sistor  connected  between  the  transmit  and  receive  transducer  terminals. 
There  is  electrical  feedthrough,  acoustic  feedthrough  through  the  cavity 
walls,  and  acoustic  feedthrough  on  modes*  other  than  the  resonant  mode 
in  the  gas  itself.  The  feedthrough  adds  vectorially  with  the  cavity  signal. 
Therefore,  in-phase  and  out  of  phase  components  of  the  composite  mea¬ 
surement  must  be  found  and  corrected  for  parallel  feedthrough. 

Information  about  the  feedthrough  is  contained  mainly  in  the  wings 
(|  co  —  oj 0 1  >  linewidth)  of  the  resonance  line.  The  signal  amplitude  in  the 
wings  depends  mainly  on  the  distance  from  the  center  frequency;  there¬ 
fore  the  resonance  signal  phase  and  amplitude  in  the  wings  can  be  calcu¬ 
lated  from  a  rough  knowledge  of  u0.  The  vector  difference  between  the 
calculated  resonance  signal  and  the  measured  composite  signal  in  the 
wings  can  then  be  interpreted  as  the  feedthrough.  This  feedthrough  can 
then  be  vectorially  subtracted  from  the  measured  composite  signal  within 
one  linewidth  of  <oQ  to  obtain  a  refined  value  of  uQ  and  linewidth  6 .  One 
can  repeat  this  procedure  to  obtain  successively  better  estimates  of  feed¬ 
through  and  resonance  signal.  Most  of  the  time  the  feedthrough  is  less 
than  5%  of  the  signal  and  one  iteration  is  sufficient  to  obtain  and  6  to 
within  the  uncertainty  of  the  measurement  (see  Fig.  2). 


*Note  that  the  resonant  mode  is  built  up  by  a  factor  of  Q  (~100)  compared 
to  a  mode  whose  center  frequency  is  far  from  the  operating  frequency. 
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The  measured  phase  at  maximum  amplitude  is  an  easily  identifi¬ 
able  point  on  the  resonance  curve.  The  phase  shift  due  to  the  resonance 
at  this  frequency  is  90  .  The  systematic  phase  shift  may  then  be  elim¬ 
inated  by  subtracting  the  phase  shift  at  ooQ  from  each  measured  phase 
and  adding  90  .  In  the  presence  of  feedthrough,  electrical  or  acoustic, 
the  maximum  amplitude  frequency  is,  in  general,  the  least  affected 
characteristic  of  the  resonance  curve. 


There  is  an  optimum  way  to  use  information  about  amplitude  to 
obtain  the  approximate  resonant  frequency  and  linewidth.  The  signal  - 
to-noise  is  best  at  maximum  amplitude.  However,  the  amplitude  is  not 
sensitive  at  wQ  changes  in  frequency.  The  amplitude  at  the  half  power 
frequencies  is  more  sensitive  to  frequency  changes.  Therefore,  two 
amplitudes  on  either  side  of  the  maximum  amplitude  are  used  to  estimate 
the  resonant  frequency  (~  maximum  amplitude  frequency). 


Working  relations  can  be  derived  once  the  line  shape  is  known. 
Asstuning  the  line  is  Lorentzian,  the  measured  amplitude  should  vary 
with  frequency  as 


where 
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The  half  width  (5  )  may  be  determined  from  points  near  the  half  power 
frequencies  of  the  line  in  question  (i.  e.  ,  0.  707  AQ  pts)  from  each  side 
of  the  center.  Putting  the  amplitudes  and  corresponding  frequencies  of 
these  two  points  in  Eq.  (18),  we  can  solve  for  5  ,  obtaining 
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where  plus  and  minus  signs  refer  to  either  side  of  the  center.  Recalling 
that  the  amplitude  is  not  sensitive  to  frequency  change  right  near  the 
center,  we  can  use  6  and  Aq  to  approximate  coq 


The  first  approximation  to  the  phase  at  maximum  frequency  is  then  the 

measured  phase  0  at  u  . 

wo  o 

The  system  phase  shift  is  in  0  wo  as  well  as  all  the  cavity  signal 
phase  shifts.  Subtracting  0WO  from  each  measured  phase  shift  eliminates 
the  system  phase  shift.  In  fact,  adding  90°  to  this  phase  difference,  the 
phase  between  the  pressure  incident  on  the  cavity  and  that  transmitted 
through  the  cavity  is  obtained.  This  resonance  phase  shift  is  calculated 
in  radians  from 
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resonance  phase  angle  interpolated  to  wq 
from  nearby  measured  phases  0  ^  and  0  ^ 
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The  systematic  phase  shift  is  not  entirely  independent  of  frequency. 
Frequency  variation  of  phase  shift  is  small  over  a  linewidth,  therefore,  it 
is  adequate  to  calculate  the  difference  between  phase  at  wQ  and  phase  at 
every  other  frequency  theoretically.  The  small  phase  variation  which  does 
rise  is  subtracted  from  the  phase  data  point  by  point. 

The  largest  mechanical  phase  shift  is  that  due  to  multiple  reflec¬ 
tion  of  the  sound  waves  between  the  foil  side  of  the  coupling  hole  and  the 
cavity.  The  phase  shift  due  to  transit  through  the  coupling  holes  is  the 
0hofEq.  (11). 
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In  general,  the  electronic  coupling  and  isolation  are  carefully 
designed  to  minimize  RC  phase  shifts.  However,  a  certain  amount  of 
bandpass  filtering  is  required  to  optimize  the  accuracy  of  the  phase 
and  amplitude  measurements.  A  tuned  amplifier  with  a  Q&  of  5  is 
used  in  the  measurements  to  date.  The  amplifier  is  set  at  the  maxi¬ 
mum  amplitude  frequency  and  the  measured  phases  at  other  frequencies 
are  corrected.  The  phase  shift  due  to  the  tuned  amplifier  is 

co 

co  _ o_ 

(0  CO 

A0t  =  tan"  — — q -  .  (21) 

a 

A  similar  correction  to  the  amplitudes  is  also  possible,  but  not 
usually  important.  Therefore  we  have  considered  only  series  phase 
shifts.  Series  phase  shifts  simply  add  to  the  desired  resonance  phase. 
Acoustical  and  electrical  feedthroughs  pass  through  channels  parallel 
to  the  desired  signal  channel.  Therefore,  the  feedthrough  must  be 
vectorially  subtracted  from  the  signal  to  obtain  the  resonance  curves. 

We  must  be  concerned  with  electrical  feedthrough  by  direct  electrical 
radiation  from  the  transmit  transducer  to  the  receive  transducer, 
acoustical  feedthrough  around  the  cavity  walls  and  mode  coupling  sig¬ 
nals  through  the  gas  path  itself.  Electrical  feedthrough  is  frequency 
independent.  However,  the  phase  of  the  acoustical  feedthrough  will  vary 

by  about  ^ir  ^  or  typically  2°  over  the  range  of  frequencies  of  interest 

for  a  given  mode. 

A  phase  diagram  showing  the  various  components  which  make 
up  the  composite  signal  (Fig.  4a)  is  helpful  in  understanding  the  method 
of  separating  series  and  parallel  signal  components.  The  lower  phase 
diagram  shows  how  the  tip  of  the  receive  pressure  vector  traces  a  circle 
counterclockwise  as  the  frequency  increases.  The  point  is  offset 
from  90°  by  the  amount  (j>  due  to  the  transit  time  in  the  coupling  holes. 

The  zero  phase  point  was  chosen  to  be  the  phase  of  the  electrical  signal 
and  the  driving  transducer.  The  electrical  feedthrough  in  these  coordinates 
has  a  phase  of  zero  degrees.  The  phase  angle  of  the  acoustic  feedthrough 
via  the  cavity  walls  can  be  determined  from  the  sound  speed  in  the  wall 
material  and  I  ^  the  effective  feedthrough  path.  The  wall  feedthrough  has 
a  phase  of  about  18  for  the  001  mode  in  the  present  cavities.  The  acoustic 
feedthrough  carried  on  nonresonant  modes  similarly  has  a  phase  somewhat 


*Phase  diagrams  from  measured  modes  are  shown  in  our  previous  report. 
Ref.  3. 
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greater  than  180°  for  frequencies  near  the  001  mode.  Both  acoustic 
feedthrough  signals  have  phase  variations  of  about  1%  over  a  linewidth. 
The  limits  of  the  phase  variations  of  the  acoustic  feedthrough  are  shown 
in  Fig.  4b  by  the  dotted  line  vectors.  In  addition,  the  acoustic  feed¬ 
through  vectors  increase  phase  as  the  drive  frequency  is  increased. 
Thus,  the  composite  feedthrough  can  be  expected  to  have  any  phase 
between  0-720°  for  lines  used  in  the  present  experiment.* 

The  typical  composite  feedthrough  vector  for  the  001  mode  is 
shown  in  Fig.  4b.  Also  shown  is  the  cavity  signal  vectors  at  u  -  «0 
and  at  the  half  power  frequencies.  All  these  vectors  are  shown  as 
dotted  lines.  In  addition,  the  reference  phase  has  been  shifted  in  ac¬ 
cordance  with  the  system  phase  shift.  The  solid  lines  show  the  mea¬ 
sured  composite  signal  which  is  the  sum  of  the  feedthrough  and  the 
signal. 


The  measured  center  frequency  is  pulled  to  frequencies  slightly 
higher  than  the  true  cavity  center  frequency.  A  feedthrough  vector  at 
30  instead  of  the  160  phase  shown  in  Fig.  4a  would  lower  the  mea¬ 
sured  center  frequency.  The  maximum  pulling  of  the  center  frequency 
occurs  for  phases  of  0  or  180  .  The  center  frequency  is  pulled  at  most 
an  amount  equal  to  the  feedthrough  amplitude  divided  by  maximum  signal 
amplitude  times  the  linewidth.  Notice  there  is  no  effect  on  center  fre¬ 
quency  for  feedthrough  vectors  with  phases  of  90  or  270  .  Similar 
remarks  hold  for  the  linewidth.  Generally,  the  feedthrough  causes  less 
than  5%  broadening  (or  narrowing)  of  the  linewidth  and  shifts  center 
frequency  by  less  than  5  Hz.  The  frequency  dependent  phase  shifts 
lead  to  ~0.  5%  effects. 

The  feedthrough  is  determined  from  a  few  measurements  (~10  6) 
far  from  the  center  frequency  and  on  either  side  of  the  line.  Only  two 
measurements  on  the  line  wings  are  required  to  obtain  feedthrough.  How¬ 
ever,  several  points  are  measured  for  a  few  lines  when  a  new  transducer 
or  cavity  is  introduced.  The  multiple  measurements  in  the  line  wings  are 
spot  checks  on  the  assumption  of  a  linear  frequency  dependence  of  feed¬ 
through,  phase  and  amplitude. 

The  line  shape  due  to  the  cavity  resonance  cannot  be  approximated 
by  a  Lorentzian  line  shape  in  the  wings  as  is  indicated  in  the  final  paragraphs 


*Most  of  the  lines  used  for  data  acquisition  in  the  present  effort  are  below 
the  third  longitudinal.  The  acoustic  feedthrough  vector,  in  general,  can 
be  n360  if  the  highest  mode  used  in  on. 
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Fig.  4a.  Composite  measured  signal  polar  diagram. 
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of  Appendix  A.  The  resonance  signal  drops  to  about  q  of  its  maximum 

value  AQ  and  then  varies  slowly  with  frequency.  The  signal  due  to  the 
cavity  resonance  and  the  feedthrough  is  usually  of  the  same  order  of 
magnitude  in  the  line  wings.  The  electronic  feedthrough  Vj  is  there¬ 
fore  the  vector  difference  between  the  measured  amplitude  vector  A 
and  the  vector  due  to  the  resonance 


V.  =  (A  cos  0  -  A  cos  0  }  +  j  (A  sin  0  -  A  sin  0  )  (22) 

I  s  s  s  s 


where  A  is  the  measured  amplitude  and  0  the  measured  phase.  The 
signal  due  to  the  cavity  resonance  Ag  and  the  associated  phase  shift  0 
are  calculated  in  the  first  iteration  from  u>0  and  6  from  Eqs.  (19a) 
and  (19b.  "Exact11  natural  line  shapes  are  given  for  our  cavity  in 
Appendix  A.  An  expression  valid  in  the  line  wings  where  |  u-w0  |  >>  6 
and  sufficiently  accurate  for  the  present  purposes  is  for  the  mode  (pn) 


(23) 


where  <»)p0  is  the  pth  radial  mode  center  frequency,  o>on  is  the  nth  longi¬ 
tudinal  mode  center  frequency,  6  the  line  width  from  Eq.  (19a).  For  the 
longitudinal  modes  <0pO  =  0  and  uon  is  obtained  from  Eq.  (19b).  For  the 
radial  modes  n  =  1  and  Wp0  is  obtained  from  Eq.  (19b)  and  measurement 
on  the  pth  radial  mode  in  question.  The  nth  radial  and  pth  longitudinal 
mode  center  frequencies  are  used  with  the  measured  5  for  compound 
modes  (pn).  The  phase  0  is  just  the  ratio  of  the  imaginary  to  real  parts 
of  the  denominator  of  Eq.  (23). 

The  feedthrough  program  calculates  A  and  0  from  Eq.  (23).  The 
true  feedthrough  is  then  obtained  by  putting  line  wing  data  and  Ag  and  0 
into  Eq.  (22).  This  feedthrough  signal  is  then  used  to  correct  the  ampli¬ 
tude  and  phase  of  the  composite  signal  measured  near  the  half  power  points. 
The  corrected  phase  versus  frequency  near  the  half  power  points  then  pro¬ 
vides  the  true  linewidth  and  center  frequency  of  the  cavity  resonance. 

^Corrected  for  systematic  phase  shifts  to  the  first  approximation. 


44 


AEDC-TR -74-33 


3.  5  SUMMARY  DATA  REDUCTION  PROCEDURE 

The  greater  part  of  Section  3  has  been  devoted  to  those  contri¬ 
butions  to  the  measured  line  shape  which  must  be  removed  to  obtain 
accurate  (0.  1%)  gas  properties  from  the  cavity  technique.  The  first 
problem  is  to  eliminate  feedthrough  the  systematic  phase  shifts.  The 
initial  phase  shift  is  set  equal  to  zero  at  the  apparent  maximum  ampli¬ 
tude  point  with  a  reference  phase  shifter.  A  rough  value  of  the  line- 
width  6  and  the  center  frequency  o>Q  is  obtained  from  Eqs.  (19a)  and 
(19b),  neglecting  feedthrough.  The  rough  values  of  6  and  w  are  then 
put  in  Eq.  (24)  to  give  the  signal  Ag  in  the  wings.  The  feedthrough  is 
then  given  by 


where  A  is  the  complex  amplitude  data  in  the  wings,  Ag  is  the  signal 
in  the  wing  and  AQ  is  the  maximum  voltage  amplitude.  R  and  1  are  the 
real  and  imaginary  parts  of  the  feedthrough. 


The  feedthrough  is  assumed  to  vary  linearly  with  frequency 

and  the  data  in  the  wings  is  used  in  Eq.  (25)  to  obtain  the  feedthrough 

constants  r  ,  r. ,  j  and  j, .  These  are 
o  1  -o  *1 


and 


r  =  r  +  r,  w 
o  1 


ijjU  . 


(25) 


(26) 


Linear  frequency  dependence  means  that  the  feedthrough  consists  of 
electrical  pickup  plus  acoustic  feedthrough  with  no  feedthrough 
resonances. 

Once  the  feedthrough  constants  are  obtained  a  better  value  of  6 
and  <*>o  are  obtained  from  the  data  within  2  6  of  u.  The  second  approxi¬ 
mation  to  and  6  is  obtained  by  fitting  the  equation 

«  =  6  (  |  a|  2  -  1)1/2  +  «o  (27) 

to  the  data.  The  quantity  a  is  then  defined  by 
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A 

a  =  —  -  r  +  ij 
o 


C 


w  -  a)  +  i6 
o 


(28) 


and  Eq.  (27)  follows  from  the  assumption  that  the  line  shape  is  Lorentzian 
within  2  6  of  u>Q.  The  variable  w  is  fit  to  a  linear  function  of  (|  a  |  ^-1)*/^. 
This  process  is  repeated  until  consistent  values  of  6  and  w  are  obtained. 


For  large  feedthrough  amplitudes. 


2  2 
R  +  I 


>  0.  3,  the  phase  shift  (<pQ) 


at  «  must  be  readjusted  after  the  first  iteration  using  the  improved  esti¬ 


mate  of  «0  and  the  0  versus  u  data, 
with  no  0  Q  adjustment  are  required. 


Most  of  the  time  two  or  three  iterations 


The  feedthrough  constants  are  then  used  to  correct  the  original 
data.  That  is,  r  is  subtracted  from  the  in-phase  signal  component  and 
j  from  the  out-of-phase  signal  component  point  by  point.  Thus  cot  0  at 
each  data  point  within  2  6  is  then  found  from  the  ratio  of  imaginary  to 
real  parts  of  the  feedthrough  corrected  signal. 


The  final  step  consists  of  making  a  fit  of  the  measured  frequencies 
to  a  linear  function  of  cot  0 


gj  =  6  cot  0  +  w  .  (29) 

'  o 

The  raw  linewidth  6  and  center  frequency  w  data  are  the  constants  of 
this  fit. 


The  coupling  hole  corrections  are  then  applied  to  the  linewidth, 
subtracting  fia  (Section  3.  2)  and  correcting  for  0  (Section  3.  1).  The 
acoustic  transmission  losses  Anp  (p.  14)  are  subtracted.  The  resulting 
linewidth  contains  only  information  on  viscosity  and  thermal  conductivity. 
This  linewidth  may  be  identified  with  the  6  -  of  Eq.  (4),  plus  the  correc¬ 
tion  factors  gyj  and  g^  from  Eq.  (14). 

The  sound  speed  is  then  obtained  from  <oQ  =  2irf0  and  Eq.  (3). 

The  sound  speed  defined  by  Eq.  (3)  is  shown  in  Appendix  A  to  be  the  real 

part  of  r— .  Since  the  internal  modes  (not  including  chemical  degrees  of 
*1 

freedom)  do  not  affect  the  isothermal  compressibility  (K.^,) 

pKT  ZRT  j^1  +  Z  (?7  jT]  (30) 
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the  real  part  of  7  is  obtained  from 


y 


R 


(31) 


Here  Z  is  the  compressibility  factor  from  pvt  measurements,  p  is  the 
density,  T  the  temperature,  and  R  the  ideal  gas  constant.  The  real  part 
of  7 ^  is  precisely  the  quantity  needed  to  determine  the  Prandtl  number 
from  the  linewidth.  However,  the  specific  heat  at  constant  pressure  is 
necessary  to  obtain  the  thermal  conductivity  from  the  Prandtl  number. 

As  long  as  the  imaginary  part  of  7"  is  small  this  specific  heat  is  given 
by 


C 

P 


(32) 


This  approximation  was  used  for  all  the  thermal  conductivity  data  of  the 
present  report. 


Specific  modes  must  be  chosen  to  show  how  the  transport  properties 
are  obtained  from  the  linewidths.  The  expressions  for  the  linewidths  for 
different  modes  are  obtained  by  putting  the  appropriate  radial  and  longi¬ 
tudinal  mode  numbers  into  Eq.  (4).  For  convenience,  the  normalized 
radial,  compound,  and  longitudinal  linewidth 


10 


1/2 


6  = 


11 


01 


1/2 


1/2 


TO 


11 


01 


respectively,  are  introduced.  One  can  eliminate  viscosity  in  favor  or 
Prandtl  number  (P^)  using  the  longitudinal  and  compound  mode  linewidths 
to  get 


5^  I1  -  V01)J  -R'  [1  +  V10)] 


7^-i 


(33) 


(1  +  R')[l  +  gjJlO)]  -  (1  +  2 R)  1/2  l1  +  \(0l)]fi 


R_ 

' 

L 


*See  Section  3.  2  for  definition  of  g  and  g^ 


47 


AE  DC-T  R-74-33 


Similar  equations  for 


TV1 

r»  1/2 


can  be  obtained  from  the  01  and  11  or  11 


and  10  modes.  The  value  of  the  coupling  hole  radius  is  varied  around 


7 

its  nominal  value  until  the  parameter  — 

P 

hole  correction  6  a  and  the  parameter  g 


is  self-consistent. 


are  functions  of 


The 


a. 


Once  a  consistent  value  of  Prandtl  number  is  obtained  the  viscosity- 
can  be  calculated  from  any  of  the  measured  linewidths.  In  particular,  the 
viscosity  can  be  calculated  from  the  longitudinal  line-width  using 


UL 

P 


1/2 


tt  r  6 


7  - 1 

R'  [l  +  gn(0l)J  +  (1  +  R')  l1  +  gX<01>] 


(34) 


Notice  that  the  Prandtl  number  does  not  depend  on  the  density. 

All  the  effects  of  density  in  Eq.  (33)  are  in  measured  quantities.  Since 
it  is  simple  to  determine  the  required  linewidth  under  identical  pressure- 


volume -temperature  conditions,  the  consistency  of 


^R-1 

■O  1/2 


was  chosen  to 


correct  the  hole  radius.  The  viscosity,  on  the  other  hand,  depends  on 
the  density.  The  density  is  determined  from  the  measured  pressure  and 
the  measured  temperature  and  available  pressure-volume-temperature, 

Refs.  26-29,  in  the  present  work.  (The  pressure  measurement  is  the 
weakest  link  of  the  present  work.  The  pressures  were  measured  to  ~0.  1%.  ) 
There  are  many  other  ways  in  which  the  self-consistency  of  redundant  cavity 
measurements  can  be  used  to  check  the  various  correction  terms.  However, 
only  the  consistency  of  the  Prandtl  number  was  used  in  the  present  work. 


The  reader  has  undoubtedly  come  to  the  conclusion  at  this  point  that 
there  is  a  lot  of  arithmetic  between  the  measurement  and  the  desired  prop¬ 
erties  with  the  cavity  technique.  However,  the  cavity  measurement  provides 
the  same  information  obtained  from  simultaneous  parallel  plate  thermal 
conductivity  measurements,^  oscillating  disk  measurements, 32  pjus 
a  Joule  Thompson^  measurement.  In  addition,  the  acoustic  cavity  can 
supply  relaxation  time  measurements.^  The  response  of  the  cavity  can 
be  calculated  in  much  greater  detail  than  the  response  of  the  above  con¬ 
ventional  measurements.  Also,  the  difficult  problem  of  measuring  small 
temperature  and  pressure  changes  is  bypassed.  Therefore,  the  working 
equations  and  the  type  of  measurements  required  are  simpler  for  the  cavity 
technique  than  for  any  other  set  of  measurements  one  might  envision  to 
measure  the  same  quantities. 
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SECTION  IV 

EXPERIMENTAL  APPARATUS  AND  RESULTS 

4.  1  SUMMARY  OF  EXPERIMENTAL  PROBLEMS  WITH 
GAS -FILLED  ACOUSTIC  CAVITIES 

There  has  been  a  great  deal  of  effort  expended  on  resonant 
acoustic  cavities  for  the  measurement  of  properties  of  gases. 35,  12,  36 
Usually,  resonators  dimensionally  large  compared  to  the  cavities  of 
the  present  work,  were  used.  The  small  dimensions  of  the  present 
cavities  were  chosen  originally  to  make  the  cavity  linewidth  large, 

~200  Hz.  This  made  the  stability  requirements  on  the  frequency  re¬ 
alistic.  For  instance,  to  measure  a  200  Hz  linewidth  to  0.  01%  directly 
requires  0.  02  Hz  frequency  determination  in  ~20  kHz.  An  additional 
reason  for  using  small  cavities  is  to  minimize  the  volume  of  material 
which  must  be  maintained  under  constant  high  temperature  and  high 
pressure  conditions. 

However,  previous  work  has  employed  large  cavities  due  to  the 
difficulty  of  eliminating  acoustic  feedthrough  via  the  solid  walls  of  the 
cavity.  Several  approaches  to  this  problem  have  been  used  by  ourselves 
and  other  workers.  Magnetically  driven  pistons, 36  flexural  ceramic 
transducers, 35  membrane  coupled  transducers  and  several  other 
transduction37  schemes  have  been  used  to  couple  more  efficiently  to 
the  gas  than  to  the  walls.  Isolation  of  the  cavity  walls  can  be  accom¬ 
plished  by  mounting  the  transducers  on  soft  members 38  such  as  O-rings 
and  minimizing  the  mechanical  contact  between  the  transducer  and  the 
cavity  walls.  Fritsche^  used  the  fact  that  flexure  of  a  thin  membrane 
provided  good  isolation  between  the  transducer  and  the  cavity  walls. 
Following  Fritsche,  we  expended  a  great  deal  of  effort  during  the  initial 
part  of  the  present  work  to  realize  a  small  membrane  isolated  trans¬ 
ducer  such  as  is  shown  in  Fig.  5a.  However,  the  small  size  of  the 
required  transducer  dimensions  meant  that  a  small  diameter  membrane 
had  to  be  used.  The  decreased  diameter  led  to  stiffening  of  the  flexure 
and  consequently  acoustic  feedthrough  and  membrane  resonances. 

Nevertheless,  our  experience  indicated  that  membrane  flexure 
was  the  most  promising  isolation  technique.  Therefore,  electrostatic 
transducers39  were  tried.  These  transducers  were  mounted  facing 
the  cavity  with  no  tension  on  the  dielectric  membrane.  The  elastic 
compliance  of  these  transducers  was  due  to  a  central  vent  in  the  trans¬ 
ducer  backup  and  the  gas  trapped  between  the  backup  and  the  transducer 
membrane.  These  transducers  are  isolated  to  a  level  below  our  detec¬ 
tion  capabilities.  They  have  a  uniform,  flat  response,  between  5  kc-250  kc. 
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That  is,  they  do  not  introduce  additional  resonances  themselves.  How¬ 
ever,  the  end  wall  impedance  is  difficult  to  control  and/or  calculate 
when  mounted  facing  the  cavity  as  shown  in  Fig.  5b.  This  problem 
was  solved  by  isolating  the  transducer  from  the  cavity  with  a  solid  end 
wall  containing  a  small  coupling  hole.  Although  the  calculation  of  the 
losses  of  energy  from  the  cavity  into  the  coupling  hole  is  involved,  it 
can  be  carried  out  with  great  precision.  In  addition,  the  end  wall  im¬ 
pedance  can  be  varied  by  simply  changing  the  size  of  the  coupling  hole. 

Finally,  the  electronic  readout  of  the  phase  and  amplitude  versus 
frequency  must  be  done  in  sufficient  detail  to  eliminate  interference  from 
feedthrough  and  thermal  drift.  Also,  problems  with  resonances  in  the 
transducer  and  cavity  walls  can  be  detected  only  from  detailed  line  shape 
measurement.  The  first  technique  used  consisted  of  setting  severed, 
frequency  points  around  a  resonance  and  measuring  the  phase  and  am¬ 
plitude  at  each  point.  The  most  serious  problem  with  this  approach  is 
that  the  frequency  must  be  stable  over  the  period  of  the  measurement, 
typically  to  0.  02  Hz  in  20  kHz  (1  pt  in  10^).  A  factor  of  three  variation 
in  center  frequency  is  necessary  to  make  a  measurement  with  the  cavity 
technique.  Therefore,  parts  in  10^  stability  is  difficult  to  achieve.  The 
proliferation  of  digital  frequency  synthesizers  during  the  past  two  years 
has  solved  this  problem. 

However,  one  is  left  with  thermal  drift  of  the  cavity  center  fre¬ 
quency  due  to  the  temperature  dependence  of  sound  speed.  The  precise 
calibration  of  phase  shift  to  0.  01%  is  also  difficult.  Several  modulation 
techniques  were  developed  in  which  two  waves  spaced  a>m  ~  6  apart  and 
centered  on  a  frequency  equal  to  the  center  frequency  of  the  mode  in 
question  were  simultaneously  sent  through  the  cavity.  Finding  the  half 
power  |450  points)  was  reduced  to  varying  wm  until  a  null  was  observed. 
The  90  phase  shift  was  automatically  assured  by  the  theory  of  the  modula 
tion  process.  However,  the  modulation  techniques  were  not  suited  to 
collecting  sufficient  information  to  process  out  interfering  signals.  There 
fore,  a  line  sweep  was  found  necessary  to  determine  feedthrough,  etc. 
Then,  additional  measurements  of  the  modulation  techniques  could  also 
be  used  to  obtain  extremely  precise  (pts  in  10®)  linewidth  determinations. 
However,  all  the  data  in  the  present  report  were  taken  with  the  line  sweep 
technique.  The  modulation  device  was  used  only  to  get  a  quadrature  cali¬ 
bration  for  the  phase -lock  amplifiers. 

4.  2  MECHANICAL  HARDWARE 

The  cavity  is  mounted  inside  a  pressure  vessel.  Most  of  the 
measurements  were  made  in  an  Aminco  pressure  cell.  The  pressure 
cell  had  two  thermocouple  ports,  two  electrical  feedthroughs,  and  a  gas 
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No  dean  modes  above  first  longitudinal  (19000  KMz)  obtained, 
(a)  MEMBRANE  ISOLATION 


METALLIC 

FILM 


AIR  GAP 


BACKUP  ELECTRODE 

PT 


BOND 


MYLAR  MEMBRANE 


CAVITY 


First  twenty  longitudinal  modes  observed- no  radial  modes. 


(b)  ELECTROSTATIC 


Fig.  5.  Techniques  for  exciting  acoustic  cavities  with  good  wall  isolation. 
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inlet  and  outlet.  A  thermistor  was  inserted  through  one  of  the  thermo¬ 
couple  ports  on  all  measurements  below  400  K,  Two  high  pressure 
gages  0-1500  and  0-7500  psia  -  Bourdon  gages  (Heise)  -  were  mounted 
on  one  of  the  gas  ports.  A  low  pressure  gage  0-100  (also  Heise)  and 
a  mercury  manometer  were  attached  through  a  second  port.  The  sec¬ 
ond  port  was  used  to  load  and  evacuate  the  high  pressure  cell.  The 
details  are  given  in  Fig.  6. 

A  second  cell  was  set  up  for  high  temperature  (700°K),  high 
pressure  (500  atm)  operation.  The  cell  material,  Hastalloy  C,  was 
chosen  specifically  for  high  pressure-high  temperature  operation.  The 
cell  volume  was  minimized  for  safety.  The  cavity  slip-fit  into  the 
pressure  cell,  as  indicated  in  Fig.  7.  The  electrical  signal  lines  were 
brought  in  through  the  1/8  inch  ID  high  pressure  tubing.  The  high  pres¬ 
sure  lines  from  the  cell  are  ~25  cm  long  and  made  of  Hastalloy  C.  The 
high  pressure  electrical  feedthroughs  are  made  at  low  temperature  on 
a  high  pressure  tee  at  the  end  of  the  Hastalloy  C  high  pressure  lines. 

The  gage,  vacuum  and  source  connections  were  common  to  the  first  cell. 

A  high  temperature  quartz  oven  was  used  for  measurements  in 
nitrogen  and  air.  In  addition,  elevated  temperature  measurements  were 
made  in  atmospheric  air  with  a  cavity  inside  a  heating  element  wrapped 
in  insulating  material.  The  quartz  oven  shown  in  Fig.  8  had  a  cavity 
mounted  over  three-quarters  of  the  way  down  a  4  inch  diameter  quartz 
tube.  The  cavity  was  mechanically  supported  by  two  steel  tubes  which 
in  turn  were  rigidly  fastened  to  the  steel  end  flange  of  the  furnace.  The 
electrical  signals  were  brought  into  the  oven  through  these  steel  tubes 
in  quartz  tubing.  The  vacuum  seal  of  the  electrical  lines  and  mechanical 
support  for  the  cavity  were  accomplished  at  the  cool  oven  end  flange. 

The  18  inch  long  quartz  tube  was  immersed  14  inches  into  a  fur¬ 
nace.  The  last  4  inches  of  the  quartz  tube  was  exposed  to  the  room 
temperature  air.  The  exposed  parts  of  the  quartz  and  the  end  flange  were 
always  sufficiently  cool  to  touch.  Four  stainless  steel  baffles  were  placed 
along  the  length  of  the  quartz  oven  to  prevent  radiation  and  convection. 
Since  the  quartz  and  steel  pieces  were  only  fastened  together  at  one  end 
no  thermal  expansion  problems  were  encountered.  However,  the  central 
steel  wire  conductor  of  the  signal  feedthroughs  was  coiled  to  avoid  thermal 
expansion  against  the  transducer  membrane. 

The  sound  speed  was  taken  as  the  true  measure  of  temperature 
in  1  atm  high  temperature  measurements.  However,  theremocouples 
were  inserted  into  a  hole  about  halfway  into  the  wall  of  the  acoustic 
cavities.  Unfortunately,  it  was  not  possible  to  measure  temperature 
gradients  in  the  walls  of  the  cavity  due  to  time  limitations.  However, 
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Fig.  7.  Hastalloy  high  pressure  cell  for  high  pressure  use. 
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Fig.  8.  Quartz  oven  for  low  pressure-high  temperature  use. 
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all  the  ovens  built  around  the  cavities  had  at  least  2  inches  of  alumina 
cotton  between  the  test  volume,  and  the  heating  elements.  In  addition, 
at  least  4  inches  of  similar  insulating  material  was  packed  on  either 
side  of  the  heaters  as  indicated  in  Figs.  1  through  8.  The  construction 
was  designed  to  minimize  gradients  and  to  stabilize  the  test  zone  against 
power  fluctuations  in  the  heater.  In  fact,  all  the  ovens  used  required 
the  better  part  of  an  8-hour  day  to  come  to  equilibrium.  Equilibrium  was 
judged  by  the  lack  of  drift  in  the  phase  null  at  resonance.  Temperature 
fluctuations  of  1  pt  in  10^  could  be  detected*  this  way. 

The  quartz  oven  could  be  evacuated  and  filled  with  any  desired 
gas.  Mainly,  nitrogen  and  air  were  investigated  in  the  quartz  oven.  Air 
at  atmospheric  pressure  was  also  investigated  in  the  oven  assembly  shown 
in  Fig.  9.  This  cavity  always  had  a  gas  leak  path  to  atmospheric  air.  In 
addition,  the  cavity  body  was  segmented  to  allow  the  coupling  hole  correc¬ 
tion  to  be  studied  for  a  given  set  of  end  plates  on  different  length  cavities. 
Discontinuities  of  up  to  0.  0001"  may  have  existed  between  the  radii  of  the 
various  segments  of  the  cavity  plates.  However,  the  data  agree  with  the 
quartz  cavity  data  to  0.  5%.  The  temperature  determinations  in  the  seg¬ 
mented  cavity  were  accomplished  in  a  way  identical  to  those  used  in  the 
quartz  ovens. 

The  acoustic  cavity  assemblies  used  during  the  present  work  were 
constructed  in  several  radial  pieces  bolted  together  as  in  Fig.  10.  The 
acoustic  cavity  itself  is  formed  by  the  central  cylindrical  piece  and  the 
end  plates.  The  gas  is  allowed  to  leak  into  the  cavity  between  the  end 
plates  and  the  side  walls.  Thus,  it  is  necessary  to  pull  a  vacuum  on  the 
cell  or  oven  containing  the  cavity  and  load  to  a  high  pressure  to  obtain  a 
pure  test  gas.  A  small  radial  scratch  along  the  side  wall  was  sufficient 
to  assure  good  exchange  of  gas.  Typically,  three  evacuation  and  load 
cycles  were  sufficient  to  assure  a  pure  gas  in  the  cavity.  Purity  of  the 
gas  was  checked  by  comparing  the  sound  speed  at  one  atmosphere  during 
two  successive  phases  of  the  load  cycle.  In  addition,  the  time  required 
for  the  gas  to  leak  into  the  cavity  and  come  to  thermal  equilibrium  with 
the  cavity  walls  could  be  determined  by  monitoring  the  change  of  sound 
speed  with  time  after  a  positive  pressure  change.  Typically,  impurity 
levels  of  5  pts  in  10^  were  observable  and  a  15-minute  equilibrium  time 
was  measured. 

The  end  plates  were  held  onto  a  thick  walled  (greater  than  1/4") 
cavity  head  by  small  screws  to  facilitate  initial  test  and  assembly  of  the 


*Absence  of  drift  implied  temperature  stability.  However,  the  presence 
of  drift  could  be  due  to  several  mechanisms. 
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Fig.  10.  Acoustic  cavity  and  transducer. 
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transducers.  Thick  walled  transducer  heads  were  found  necessary  to 
reduce  electrical  pickup  between  the  receive  and  transmit  transducers. 
In  addition,  such  construction  minimized  the  possibility  of  any  low 
frequency  acoustic  resonances  in  the  head  itself. 


The  end  plates  were  made  of  the  same  material  as  the  cavity 
body.  The  end  plates  were  between  0.  330  and  0.  139  cm  thick.  Thinner 
end  walls  than  ~0.  1  cm  had  plate  resonances. 


The  area  of  the  end  plate  was  ~2000  times  the  area  of  the  coupl¬ 
ing  holes.  Thus,  plate  resonances  were  extremely  efficient  in  exciting 
sound  in  the  cavity.  Any  sound  introduced  this  way  into  the  gas  had  plate 
resonances  (Ref.  40,  p.  213)  and  phase  shifts  superimposed  on  cavity 
resonances.  End  plate  resonances  become  more  troublesome  at  elevated 
temperatures  because  several  factors  conspire  to  drop  the  plate  resonance 
frequencies  as  temperature  increases.  Large  radius  cavities  (large  end 
plate  radius)  must  be  used  to  maintain  a  sufficiently  narrow  linewidth  at 
elevated  temperatures.  Also,  as  temperature  increases  the  elastic 
modulus  of  the  steel  end  plate  drops.  Thus,  lowest  frequency  end  plate 
resonance  was  designed  to  be  ~80  kHz  at  room  temperature  (~1  cm 
radius  cavity).  A  thick  (~0.  330  cm)  plate  was  used  with  a  thin  (~0.  139) 
plate  to  avoid  simultaneous  resonances  in  the  two  end  walls. 


Cavities  with  radii  ranging  from  1.  26  cm  down  to  0.  25  cm  and 
radius  to  diameter  ratios  ranging  between  0.  1  and  1  were  successfully 
employed  in  the  present  work.  However,  only  five  cavities  were  used 
to  take  the  data  presented  below.  The  material  and  dimensions  of  these 
cavities  are  presented  in  Table  4.  The  thermal  expansion  of  the  cavity 
material  is  also  shown  in  Table  4.  The  change  in  cavity  dimensions  due 
to  thermal  expansion  between  ~300  and  1000  K  amounts  to  ~1%.  Mea¬ 
surement  of  the  radial  mode  resonant  frequency  versus  temperature  in  a 
gas  such  as  argon  at  low  density  ('"T  amagat)  would  suffice  to  determine 
cavity  dimensions  versus  temperature  to  parts  in  10^.  However,  available 
thermal  expansion  data^*  was  used  to  correct  viscosity  data.  Recall  in 
Section  3.  5  that  only  the  cavity  dimensions  enter  the  Prandtl  number 


determination.  Thus,  the  Prandtl  number  measurement  (actually 
is  not  dependent  on  thermal  expansion. 


Y-1 


1/2 


The  cavity  dimensions  and  materials  of  the  present  study  were 
chosen  for  flexibility  in  the  face  of  a  program  calling  for  measurements 
over  a  wide  range  of  pressures,  temperatures  and  gases.  Future  work 
should  employ  cavities  optimized  for  more  specific  areas  of  the  pressure- 
temperature  range.  High  acoustic  impedance  materials  such  as  tungsten 


59 


Table  4.  Material  and  dimensional  characteristics  of  cavities  used  to  measure 


data  reported  below. 


Cavity  and 

Test  Volume 

Radius 

cm 

R' 

Material 

Thermal 

Expansion 

Coefficient 

Precision  (Aminco 
high  pressure  cell) 

0.  72896 

0.  81223 

5200  steel 

8.  4  x  10"6  (0-400°K) 

High  pressure  (Aminco 
high  pressure  cell) 

0.  72433 

0.  33621 

302  Stainless 
steel 

9.  6  x  10“6  (0-400°K) 

High  temperature 
(quartz  oven) 

1. 2584 

0.  81215 

II 

9-  6  x  10"6  (0-400°K) 

Segmented  (open  oven 
atmospheric  air) 

0.  62808 

0.  78039 

(varies  with  run) 

Cold  rolled 
steel  1018 

8.  4  x  10~6  (0-900°K) 

High  temperature, 
high  pressure 
(Hastalloy  pressure 
cell) 

0.  64255 

0.  37683 

. 

Invar 

0.7  x  10"6  (300°K) 

6.4  x  10"6  (600°K) 
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will  reduce  acoustic  transmission  losses  by  a  factor  of  5  or  more  for 
high  density  (greater  than  100  amagats)  measurements.  Oxidation 
resistant  nickel  base  superalloys  can  be  used  effectively  for  measure¬ 
ments  in  air  and  oxygen  at  temperatures  up  to  ~1000  K.  Finally,  the 
advent  of  digital  frequency  synthesizers  at  low  frequencies  together 
with  tracking  phase-lock  amplifiers  means  that  large  cavities  (~3") 
can  be  used  without  compromising  electronic  stability  requirements. 
Larger  cavities  simplify  the  hole  correction  problem  and  allow  high 
cavity  Q  to  be  obtained.  Large  cavities  have  advantages  for  low  pres¬ 
sures  (under  100  atm)  and  high  temperature  measurements.  Neverthe¬ 
less,  any  special  optimization  introduces  problems.  High  pc  and 
oxidation  resistant  materials  are  difficult  to  procure  and  fabricate. 

Large  cavities  make  thermal  stability  and  uniformity  problems  more 
acute.  In  addition,  large  cavities  pose  difficult  pressure  cell  design 
problems  especially  above  ~200  C.  These  are  the  basic  considerations 
behind  the  mechanical  design  and  the  cavities  for  the  present  work. 

4.  3  TRANSDUCERS 

The  electrostatic  transducers,  shown  in  Fig.  10,  consist  of  a 
flexible  dielectric  between  a  flexible  conducting  layer  and  the  back  side 
of  the  cavity  end  wall.  A  dc  voltage  between  the  cavity  end  plate  (ground) 
and  the  flexible  layer  cause  a  compressive  electrostatic  force  between 
the  two  conducting  parts.  An  alternating  signal  superimposed  on  the  dc 
field  causes  the  electrostatic  force  to  fluctuate.  The  dielectric  then  also 
vibrates,  driving  the  trapped  gas  back  and  forth  in  the  coupling  hole. 
Similarly,  a  pressure  wave  in  the  coupling  hole  causes  the  dielectric 
film  to  vibrate.  The  resulting  variations  in  capacitance  of  the  electro¬ 
static  transducer  causes  an  rf  voltage  to  appear  across  the  receive 
transducer. 

In  fact,  the  operation  of  the  electrostatic  transducers  is  consider¬ 
ably  complicated  by  the  formation  of  polarization^  in  the  dielectric.  The 
dc  bias  gradually  induces  dipole  states  both  on  the  surface  of  the  dielec¬ 
tric  and  in  the  body  of  the  dielectric  layer.  The  polarization  field  opposes 
the  bias  field.  As  a  result  the  signal  of  the  new  transducer  falls  off  as  a 
function  of  time.  Different  dielectrics  show  varying  ability  to  polarize. 
Mylar  and  mica  show  100%  polarizability  so  that  the  signal  would  drop  off 
to  zero  eventually.  The  teflon  and  quartz  dielectrics  show  very  little  sig¬ 
nal  loss  due  to  polarization.  In  addition,  above  a  certain  temperature, 
which  we  have  called  the  depolarization,  no  permanent  polarization  is 
observed. 

The  time  dependent  polarization  is  a  major  problem  just  below 
the  depolarization  temperature  because  the  amplitude  change  due  to 
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depolarization  is  superimposed  on  changes  due  to  thermal  drift  as  well 
as  those  due  to  the  resonance.  There  is  also  drift  in  tile  system  phase 
shift  since  the  space  between  the  conducting  plates  increases,  as  the 
polarization  field  cancels  the  dc  bias  field,  decreasing  the  transducer 
capacitance.  Therefore,  it  is  necessary  to  operate  above  or  well  below 
the  depolarization  temperature.  A  drift  of  ~1%  in  amplitude  per  hour 
was  found  to  be  tolerable. 

Several  materials  were  tried  for  the  flexible  parts  of  the  electro¬ 
static  transducer.  The  results  are  summarized  in  Table  5.  Our  early 
work  used  mylar  films  as  the  dielectric.  The  mylar  was  available  with 
an  aluminum  coating.  A  mylar  transducer  built  up  a  polarization  field 
at  room  temperature  in  a  matter  of  minutes  on  the  receive  side  after  an 
~100  dc  bias  was  turned  on.  The  receive  side  could  then  be  operated 
without  bias  for  several  days  before  noticeable  depolarization  occurred. 
However,  the  drive  signal  (100  at  ~20  kHz)  was  sufficient  to  keep  the 
transmit  side  depolarized. 

The  aluminized  mylar  is  very  satisfactory  until  a  slight  tempera¬ 
ture  range  is  introduced.  Above  50  K,  in  the  absence  of  the  dc  bias,  the 
mylar  becomes  fully  depolarized  in  a  matter  of  minutes.  It  is  then  not 
possible  to  obtain  a  stable  signal. 

Quartz  and  mica  dielectric  membranes  are  useful  at  high  tempera¬ 
tures.  The  quartz  has  the  advantage  of  very  little  polarization.  Quartz 
films  can  easily  be  obtained  by  blowing  bubbles  in  quartz  tubing.  However, 
it  is  difficult  to  obtain  flat  quartz  films.  In  addition,  the  films  fail  by 
fracture  to  easily  be  practical.  However,  more  recent  work  which  we 
have  not  had  time  to  complete  indicates  that  much  thicker  films,  ~0.  002, 
might  be  satisfactory  from  an  acoustic  point  of  view.  Small' diameter  flat 
plates  of  such  material  are  available  from  Corning  Glass  Works. 

The  mass  of  the  movable  films  increases  the  amount  of  coupling  of 
sound  between  the  end  plate  and  the  transducer.  An  unbonded  metal  foil 
is  equally  effective  for  the  flexible  electrode  as  a  bonded  foil  or  coating. 

In  addition,  transducers  constructed  with  metal  foils  do  not  fail  by  electrode 
wear  as  quickly  as  coated  transducers.  In  addition,  unbonded  foils  have 
obvious  advantages  from  the  point  of  view  of  thermal  expansion  on  delicate 
dielectrics.  Nevertheless,  tests  on  the  large  diameter,  high  temperature 
cavity,  showed  a  great  deal  of  end  plate  excitation.  For  this  reason  the 
back  electrode  was  made  with  a  much  thinner  layer  of  vapor  deposited 
nickel -chrome. 

Mica  is  very  useful  for  high  temperature  dielectric  films.  Com¬ 
mercial  plates  ~0.  005"  thick  can  easily  be  split  into  ~5  parts  under  water. 
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Table  5.  Dielectric  and  flexible  electrode  parameters. 


Dielectric 

Electr 

ode 

Percent 

Polarizability 

Depolariz  ation 
Temperature 

Maximum 

Operation 

Temperature 

Comments 

Material 

Thickness 

inch 

Material 

Thickness 

inch 

3 

- 

250°C 

Low  activity, 
low  spurious  ‘ 
resonances 

Teflon 

0.  001 

A1  foil 

0. 0005 

Mylar 

0. 0002- 
0.  0005 

A1  coating 
Cu  foil 

io"5 

.  0003 

100 

<50 

200°C 

Depolarized  by 
100VAC  on  trans¬ 
mit  side.  Good 
activity.  Low 
spurious  reso¬ 
nances  with  A1 
coating. 

Quartz 

0.  0001- 
0.  001 

A1  foil 

Cu  foil 

Be  foil 

0.  0005 

0.  0003 

0.  0005 

1 

- 

i300°C 

Very  difficult  to 
handle.  Low 
activity. 

Mica 

0.  0005- 
0.  001 

Mo  foil 
Be  foil 
Chromel-B 

0.  0003 

0. 0005 
10'5 

100 

280 

1000°C 

Not  depolarized 
on  transmit  side. 
High  activity 
low  spurious 
resonances. 

Polarization 
greatly  enhanced 
by  H20,  C02. 
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These  films  cut  into  cm  circular  samples  and  coated  with  nickel- 
chrome  (over  a  ~0.  8  cm  diameter)  are  the  best  high  temperature 
transducers  found  to  date.  The  activity  of  the  mica  dielectric  trans¬ 
ducers  is  much  (x  5)  greater  than  transducers  with  other  dielectrics. 

This  is  possibly  due  to  the  stiffness  of  the  mica  dielectrics  causing 
more  gas  motion  around  the  coupling  hole. 

However,  the  mica  coated  with  nickel- chrome  transducer  is 
not  completely  satisfactory.  Even  after  careful  surface  preparation 
the  nickel- chrome  coating  tends  to  flake  off.  Although  the  mica 
dielectric  has  never  failed  the  transducers  do  not  survive  thermal 
cycling  above  600  K.  In  addition,  the  transducer  coating  lasts  only 
about  8  hours  at  1000  K.  The  mica  transducers  show  strong  polar¬ 
ization  effects.  The  polarization  rate  seems  to  vary  somewhat  from 
sample  to  sample  and  is  greatly  enhanced  by  H^O  and  CC^.  In  dry 
N£»  Ar  and  air  the  mica  dielectric  transducers  work  well  except  for 
a  small  temperature  range  between  200  and  280  K  where  polarization 
changes  very  quickly. 

One  of  the  difficult  problems  in  realizing  the  electrostatic 
transducers  is  obtaining  small  coupling  holes  with  diameters  of 
<0.  010".  The  hole  has  to  be  in  a  thick  enough  plate  so  that  reson¬ 
ances  in  the  plate  vibrations  are  all  well  above  the  frequencies  of 
interest.  This  means  a  0.  010"  hole  in  a  0. 100"  plate.  After  trying 
several  potential  sources  of  such  holes  we  found  laser  drilling  and 
jewel  bearing  firms  which  could  provide  apprbpriate  coupling  holes. 

A  company  called  Laser  Nucleonics  in  Waltham  has  put  a  0.  004"  hole 
in  a  0.  125"  plate  and  can  drill  smaller  holes  if  desired.  A  second 
approach  has  involved  using  orifice  jewels  (made  by  watch  jewel  manu¬ 
facturers  in  Waltham).  We  have  tried  jewels  with  holes  ranging  from 
0.  008"  to  0.  002"  in  diameter.  The  geometry  of  these  holes  is  extremely 
well  known  which  makes  orifice  jewels  particularly  attractive  for  check¬ 
ing  coupling  hole  corrections.  However,  most  of  the  coupling  holes  used 
for  data  taking  were  smooth  bore  straight  through  holes  drilled  by  con¬ 
ventional*  means,  with  diameters  between  0.  008"  and  0.  013". 

The  laser  drilled  holes  are  conical,  tapering  about  0.  002"  per 
0.  020"  of  length.  The  walls  of  the  hole  are  quite  rough  and  the  modes 
on  the  coupling  hole  cannot  be  calculated.  However,  the  laser  drilled 
holes  give  a  large  amplitude  of  cavity  excitation  with  a  small  coupling 
correction.  The  laser  drilled  holes  must  be  calibrated  to  obtain  1% 
corrections,  which  is  equivalent  to  0.  01%  in  linewidth. 


*Electric  discharge  and  cut  drill  stocks. 
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The  orifice  jewel  coupling  holes  are  smooth  bore  cylindrical 
holos.  The  jewel  is  a  piece  of  AI2O3  about  0.  030"  long  and  0.  060" 
wide  with  an  ~0.  005"  hole  drilled  in  the  center.  The  jewel  is  press 
fit  into  a  0.  060"  hole  drilled  in  the  end  plate  of  the  cavity.  The  jewel 
holes'  geometry  is  well  known  because  an  optical  comparator  photo¬ 
graph  is  supplied  with  each  jewel.  A  photograph  of  the  jewel  cross 
section  is  useful  because  the  jewel  is  transparent.  The  hole  and  the 
taper  in  the  back  end  are  clearly  visible  in  these  pictures.  Since  the 
outside  radius  is  accurately  known,  the  hole  length  and  radius  can  be 
determined  from  the  photographs  to  better  than  1%.  Thus,  the  jewels 
are  ideal  for  calibration  purposes.  It  should  be  noted  however  that 
the  output  from  jewel  coupling  holes  is  sharply  reduced  compared  to 
drilled  holes.  This  can  be  rectified  by  using  a  metallic  backing  on  the 
jewel. 


The  coupling  holes  used  for  the  bulk  of  the  data  were  uniform 
radius  (to  ~1%)  through  the  entire  plate  thickness.  The  holes  were 
cleaned  with  steel  cleaning  rods  which  can  be  obtained  from  watch 
supply  outlets.  Otherwise,  the  holes  were  used  in  the  as  drilled  con¬ 
dition.  Due  to  the  length  of  the  required  coupling  holes  0.  008"  was  the 
lowest  diameter  drilled  hole  we  could  obtain.  Much  smaller  holes  can 
be  obtained  by  slip  fitting  steel  capillary  tubes  to  larger  holes  in  the  end 
plate.  However,  several  acoustic  considerations  put  a  lower  limit  on 
the  permissible  hole  size.  Considerable  reduction  of  the  hole  correc¬ 
tion  can  be  effected  once  the  operation  of  coupling  holes  smaller  than 
0.  008"  is  studied  experimentally. 

Some  of  the  factors  limiting  hole  size  variation  can  be  seen  from 
consideration  of  signal  levels.  The  output  voltages  of  the  receive  trans¬ 
ducer  is  between  0.  2-20  p  V  in  all  the  measurements  reported  in  the 
present  work.  Transducer  sensitivity  is  about  1  |j.V/|xbar.  Therefore, 
the  pressure  wave  in  the  receive  coupling  hole  is  2  x  10“^  to  2  x  10"^ 
bars.  The  reflection  of  the  pressure  wave  in  the  cavity  approximately 
doubles  the  incident  pressure  amplitude  in  the  coupling  hole.  This  means 
that  pressure  levels  in  the  cavity  are  10-^-10-^  bars.  Due  to  the  Q  of  the 
cavity,  the  field  in  the  cavity  is  Q  times  the  exciting  field  from  the  trans¬ 
mitter  coupling  holes.  However,  the  pressure  incident  on  the  cavity  from 

the  coupling  hole  is  ~  j  times  as  large  as  the  pressure  transmitted 

into  the  cavity.  Therefore,  the  pressure  field  in  the  coupling  holes  is 
approximately  10"^- 10"*  bars!  Thus,  because  of  the  large  area  ratio 
between  the  coupling  holes  and  the  cavity  one  can  easily  get  nonlinear 
flow  in  the  coupling  holes.  Nonlinearities  in  the  coupling  hole  render 
the  hole  correction  intractable  and  cause  several  interfering  effects 
(such  as  odd  harmonics). 
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The  above  discussion  applies  for  0.  010"  diameter  coupling 
holes.  The  receive  voltage  level  of  0.  2  [iV  is  the  minimum  level 
that  one  can  use  without  intolerable  feedthrough  and  receiver  noise. 


Since  the  receive  signal  varies  with 
tude  in  the  transmitter  coupling  hole  must  be  raised  by 


■ft)' 


the  sound  pressure  ampli- 
2 


ft)' 


when 


the  coupling  hole  size  is  degreased.  It  appears  that  coupling  holes 
~0.  001"  diameter  should  be  possible  with  a  drive  level  of  10-^  and  a 
receive  voltage  level  of  0.  2  p.  V  from  the  above  considerations.  How¬ 
ever,  successful  operation  below  ~0.  005"  may  not  be  possible  due  to 
inefficiency  of  the  electrostatic  transducer  operating  into  small  holes 
and  viscous  losses  in  the  coupling  holes.  In  addition,  signal  levels 
higher  than  0.  2  |jl  V  are  necessary  for  linewidth  measurements  better 
than  ~0.  1%. 


4.4  ELECTRONICS 


Two  electronic  systems  were  used  to  drive  cavity  and  process 
the  output.  The  initial  work  in  argon  and  nitrogen  was  completed 
using  a  PAR  HR28  phase-lock  amplifier.  The  PAR  has  an  oscillator 
which  provided  the  (1  V)  drive  amplifier  signal.  The  drive  amplifier 
in  turn  provided  up  to  270  volts  rf  to  the  transmit  transducer.  The 
signal  from  the  output  of  the  cavity  was  fed  in  turn  into  the  signal  end 
of  the  PAR  instrument. 


The  signal  end  of  the  HR28  also  incorporates  an  amplifier  which 
is  gange  tuned  with  the  oscillator.  The  tuned  amplifier  has  an  adjust¬ 
able  Q  which  was  set  at  5.  The  amplifier  feeds  a  phase  detector.  The 
phase  of  signal  from  the  amplifier  is  compared  to  the  phase  of  a  refer¬ 
ence  signal  from  the  HR-8  oscillator.  A  phase  shifter  in  the  reference 
channel  allows  the  phase  detector  to  be  nulled.  The  phase  reading  on  the 
phase  shifter  is  then  the  90  complement  of  the  phase  shift  between  the 
oscillator  driving  the  cavity  and  the  signal  through  the  cavity. 

The  measured  quantity  is  the  phase  shift  between  the  reference 
voltage  inside  the  phase-lock  amplifier  and  the  voltage  at  the  input  to 
the  HR-8  phase  detector.  Examination  of  the  overall  experimental 
diagram  in  Fig.  11  indicates  that  there  are  several  electrical  elements 
between  the  pressure  wave  in  the  cavity  and  the  electrical  signals  we 
measure.  Therefore,  tests  were  carried  out  to  prove  that  these  elec¬ 
trical  elements  do  not  introduce  additional  phase  shifts  which  influence 
the  linewidth  measurement. 
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The  electrostatic  transducer  requires  a  ~100  volt  ac  signal 
and  a  variable  dc  bias  of  up  to  200  dc.  Since  the  signal  channel  tuning 
must  be  exactly  set  on  the  reference  frequency,  we  obtained  the  trans¬ 
ducer  driving  signal  by  amplifying  the  reference  signal  from  the  phase- 
lock  amplifier.  It  is  relatively  simple  to  build  an  amplifier  which  has 
a  roughly  linear  frequency  dependent  phase  shift.  One  such  device 
built  around  a  12  AU7  is  indicated  in  Fig.  11. 

The  output  of  the  drive  amplifier  might  seem  a  better  choice  for 
a  reference  than  the  HR28  internal  reference.  This  would  avoid  any 
possible  phase  shifts  in  the  amplifier.  However,  it  is  much  easier  to 
obtain  a  distortion  free  reference  voltage  at  low  voltages.  The  cavity 
and  the  tuned  amplifier  of  the  HR28  filter  out  distortion  introduced  by 
the  drive  amplifier  in  the  signal  channel. 

Most  of  the  data  was  taken  with  an  improved  electronic  system. 
The  block  diagram  for  this  system  is  shown  in  Fig.  12.  The  heart  of 
this  system  is  the  Model  129  PAR  tracking  phase-lock  amplifier.  This 
amplifier  puts  out  a  dc  voltage  proportional  to  the  in-phase  and  90  out 
of  phase  components  of  the  signal  channel  with  respect  to  the  reference 
channel.  These  two  components  of  the  signal  are  read  out  on  a  digital 
voltmeter. 

The  frequency  is  set  to  0.  001  Hz,  using  a  Rockland  digital  fre¬ 
quency  synthesizer.  The  synthesizer  frequency  must  be  checked  on  a 
counter  to  be  sure  it  is  working  properly.  The  sine  wave  is  fed  into  a 
drive  amplifier  which  provides  the  large  amplitude  required  to  drive 
the  cavity. 

A  cathode  follower  stage  was  added  to  the  drive  amplifier  to 
improve  the  high  frequency  output.  The  harmonic  distortion  of  the  drive 
signal  is  also  considerably  reduced  by  the  cathode  follower  stage.  In 
addition,  a  variable  transducer  bias  control  was  incorporated  in  the  out¬ 
put  of  the  drive  amplifier.  The  variable  bias  supply  is  required  because 
the  high  voltage  rf  is  not  sufficient  to  depolarize  the  mica  dielectric. 
Thus,  periodically  the  drive  transducer  bias  must  be  adjusted  between 
~10  and  40  volts.  The  improved  version  of  the  driver  amplifier  is  shown 
in  Fig.  13. 


The  frequency  synthesizer  and  tracking  phase-lock  amplifier 
reduce  the  time  required  to  make  a  complete  measurement  from  three 
hours  to  one-half  hour  over  the  system  in  Fig.  11.  A  three -hour  mea¬ 
surement  time  means  that  two  data  points  per  day  are  an  absolute  maxi¬ 
mum  data  rate.  Any  equipment  failure  or  temperature  instability  means 
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only  one  data  point  can  be  taken.  The  new  equipment  reduces  the  mea¬ 
surement  time  to  a  half  hour  so  that  any  problems  can  be  corrected, 
the  system  equilibrated,  with  enough  time  remaining  in  a  half  day's 
work  to  assure  at  least  one  measurement.  Four  data  points  a  day  can 
be  taken  at  one  temperature  and  different  pressures  with  the  Rockland 
frequency  synthesizer  and  PAR  129  phase-lock  amplifier. 

The  synthesizer  and  tracking  phase  meter  also  improve  the 
immunity  of  the  technique  to  line  voltage,  temperature,  pressure 
and  oscillator  frequency  drift.  During  a  measurement  of  phase  and 
amplitude  versus  frequency  for  a  single  mode  the  maximum  drift 
consistent  with  0.  1%  transport  property  measurement  accuracy  is 
0.  1  Hz,  6  x  10" ^  C,  and  0.  2  psi  for  frequency,  temperature  and  pres¬ 
sure,  respectively.  Measurement  of  a  single  mode  with  the  system  in 
Fig.  11  requires  these  tolerances  be  held  for  an  hour.  The  setup  in 
Fig.  12  requires  these  tolerances  be  held  for  only  ten  minutes.  Finally, 
the  synthesizer  allows  the  frequency  to  be  changed  by  exact  uniform 
steps  (for  instance,  10  Hz  to  0.  001  Hz).  Thus,  a  counter  is  not 
necessary.  In  addition,  the  uniform  spacing  of  frequency  points  de¬ 
creases  data  reduction  time  by  at  least  a  factor  of  three. 

The  Model  129  phase-lock  amplifier  automatically  tracks  the 
synthesizer  frequency,  eliminating  the  time  and  error  involved  in 
tuning  the  receive  amplifier  to  transmit  signal  frequency.  However, 
the  tuned  amplifier  of  the  HR28  is  still  used  and  is  set  to  the  line  center 
before  a  line  sweep  is  taken.  In  addition,  both  phase  and  amplitude 
measurements  are  done  in  the  phase-lock  mode.  This  means  that  the 
noise  immunity  is  improved.  Finally,  phase  and  amplitude  measure¬ 
ments  are  obtained  simultaneously  at  one  frequency  which  further 
reduces  the  data  reduction  time. 

This  is  especially  true  of  time-consuming  feedthrough  correc¬ 
tions.  The  0  (in-phase)  and  90  (out  of  phase)  components  of  both  the 
feedthrough  and  cavity  signal  had  to  be  calculated  from  the  phase  and 
amplitude  information  with  the  system  shown  in  Fig.  11.  The  appro¬ 
priate  feedthrough  component  could  then  be  subtracted  from  the  corres¬ 
ponding  signal  component.  The  PAR  Model  129  directly  measures  the 
0  and  90°  components. 

Although  the  line  shape  techniques  described  above  provide  a 
great  deal  of  information,  they  have  three  major  disadvantages.  They 
are:  (1)  the  measurement  takes  a  long  time,  (2)  the  temperature  must 
be  held  constant  during  this  time,  and  (3)  the  frequency  must  be  mea¬ 
sured  very  accurately  (typically  5  pt  in  10^  for  0.  01%  in  linewidth).  The 
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frequencies  and  corresponding  to  the  £  45°  phase  measurements 
are  subtracted  to  obtain  the  line  width  in  the  line  shape  technique.  Since 
the  linewidth  is  about  1%  of  f^  and  accurate  linewidths  require  much 
more  accurate  values  of  fj  and  f2-  Also,  the  temperature  must  be  held 
constant  (to  0.  01%)  for  roughly  ten  minutes  while  the  line  shape  is  being 
measured. 

Transport  property  measurements  can  in  principle  be  made  much 
faster,  and  more  accurately,  using  a  double  side -band  technique  which 
allows  center  frequency  and  linewidth  to  be  set  and  measured  separately. 

In  the  double  side-band  system  illustrated  in  Fig.  14  the  center  frequency 
(o>0)  is  sent  into  one  leg  of  a  balanced  modulator  and  a  frequency  corres¬ 
ponding  roughly  to  one-half  of  the  half  width  into  the  second  leg  of  the 
balanced  modulator.  The  output  of  the  modulator  then  consists  of  two 
equal  amplitude  sine  waves,  one  at  u>Q  +  a>m  and  the  other  at  -  o>m. 
These  two  frequencies  are  sent  simultaneously  through  the  cavity. 

The  output  of  the  cavity,  shown  also  in  Fig.  14,  shows  the  usual 
beat  pattern  due  to  the  interference  of  the  (t*>0  +.  wm)  signals.  The  center 
frequency  is  set  holding  u>m  constant  and  adjusting  wQ  until  the  envelope 
cross-overs  as  viewed  on  an  oscilloscope  are  sharp.  cjq  can  be  set  to  1% 
of  the  linewidth  this  way  in  a  matter  of  seconds.  More  accurate  line  center 
frequency  may  be  obtained  by  numerical  processing  of  wQ  versus  wm  which 
corrects  for  non-Lorentzian  line  shape,  and/or  electrical  feedthrough,  etc. 
The  line  center  so  defined  is  center  of  symmetry  of  amplitude  and  exact 
expressions  for  this  quantity  in  terms  of  cavity  geometry,  sound  speed  and 
linewidth  are  available  from  the  cavity  theory.  However,  for  0.  1%  trans¬ 
port  property  determinations  the  1%  of  linewidth  wQ  setting  and  the  assump¬ 
tion  of  a  Lorentzian  line  shape  are  adequate. 

The  setting  um  to  exactly  one-half  the  linewidth  requires  the  re¬ 
maining  electronic  elements  in  Fig.  14.  The  output  of  the  cavity  is  sent 
into  a  logarithmic  detector  which  converts  the  envelope  at  u>m  into  a 
cosine  wave  at  2  «m.  In  addition,  the  phase  of  the  detector  output  is 
twice  the  difference  between  the  cavity  phase  shift  and  <^>  _  at  u>Q  +  o>m 
and  wQ  -  wm,  respectively.  When  wm  is  one-half  the  linewidth  o>0  +  t*>m 
and  <*>0  -  u>m  are  brought  90°  out  of  phase  by  their  transit  through  the 
cavity.  The  output  of  the  logarithmic  detector  is  therefore  180°  out  of 
phase  with  the  original  «m. 

The  final  step  in  the  logic  of  setting  t*>m  equal  to  one-half  the  line- 
width  involves  generating  a  cosine  wave  at  2  a>m  but  90  out  of  phase  with 
the  original  wm  driver.  This  is  done  to  one  part  in  10^  (in-phase)  with  a 
digital  frequency  and  phase  synthesizer  described  below.  The  reference 
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Fig.  14.  Block  diagram  of  automated  linewidth  measurement  system  using  a  double  side  band  technique. 
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signal  cos  (2  wm  +  90°)  is  fed  into  the  reference  leg  of  a  balanced 
modulator  and  the  cos  [  2  tom  +  2  ($.|_  -  <£_)]  from  the  logarithmic 
detector  is  fed  into  the  other  leg.  com  is  varied  until  a  null  in  the 
balanced  modulator  output  is  obtained.  This  indicates  that  2  (<p  +  -  <p  _) 
is  180  ,  or  equivalently,  tom  is  one-half  the  linewidth. 

The  above  procedure  requires  some  time  to  describe  but  the 
setting  and  measurement  of  co0  and  u>m  are  accomplished  in  seconds. 

In  addition,  because  the  phase  is  generated  digitally  to  one  part  in 
105  the  linewidth  determinations  to  one  part  in  10^  are  possible!  Thus, 
with  a  very  simple  electronic  apparatus  this  technique  holds  promise 
for  unheard  of  precision  in  thermodynamic  and  transport  property 
measurements. 

The  key  factor  which  allows  such  precise  linewidth  measure¬ 
ment  to  be  performed  is  the  development  of  a  digital  frequency 
synthesizer  and  digitally  phase  shifted  output. 

The  problem  was  to  obtain  a  signal  source  (u>m)  which  would 
continuously  time  the  range  from  0-10  kHz  while  providing  an  ex¬ 
tremely  precise  signal  at  2  +  90°  which  would  contain  essentially 

no  even  harmonic  distortion.  These  conditions  are  met  by  the  use 
of  a  specially  modified  hybrid  circuit  manufactured  for  military  use 
by  Micro  Networks  Corporation  of  Worcester,  Massachusetts  (see 
Fig.  15).  This  hybrid  circuit  consists  of  a  seven  bit  counter,  a  1024 
bit  read-only  memory  (programmed  with  the  sine  table),  a  seven  bit 
digital-to- analog  voltage  converter  and  the  associated  reference 
voltage  supplies  and  amplifiers. 

In  operation  the  counter  is  driven  by  a  clock  signal  at  a  fre¬ 
quency  of  128  «m.  Each  pulse  of  the  clock  selects  in  sequence  one 
word  (8  bits)  of  the  memory.  The  word  is  a  number  corresponding 
to  the  amplitude  of  the  sine  wave  at  that  time  and  is  converted  to  a 
voltage  by  the  D/A  converter.  The  output  of  the  unit  is  then  a  highly 
accurate  sine  wave  whose  frequency  is  derived  from  that  of  the  clock 

_  clock  frequency 
Wm  128 

The  modifications  which  were  made  for  us  were  to  bring  out  the 
leads  from  the  5th  and  6th  bits  of  the  counter.  The  voltage  at  the  5th 
bit  changes  its  state  every  45  (that  is,  every  eighth  clock  pulse)  while 
the  voltage  from  the  6th  bit  changes  state  every  90  (every  fourth  clock 
pulse). 
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The  sine  wave  output  «m  and  the  square  waves  at  2  «m  (bit  6), 

4  u>m  (bit  5)  and  2  +  90  (exclusive  OR),  are  shown  in  the  timing 

diagram,  Fig.  1.6.  The  +  90°  square  wave  is  obtained,  for  instance, 
by  decoding  the  2  tom  and  4  um  in  an  "exclusive"  or  logic  gate.  This 
device  constitutes  the  blocks  labeled  "sine  wave"  and  "decode"  in 
Fig.  14. 


The  digital  phase  synthesizer  proved  valuable  for  calibrating  the 
quadrature  of  the  phase-lock  amplifiers  used  in  the  line  sweep  measure¬ 
ments.  The  2  «m  outputs  from  bit  6  and  the  exclusive  OR  are  90°  out  of 
phase  to  a  precision  limited  by  the  jitter  of  the  digital  logic.  This  jitter 
is  some  fraction  of  the  switching  time  (~5  nsec)  of  the  logic.  Thus,  at 
the  maximum  frequency  (2  tom  =  20  kHz)  the  device  will  put  out  the  90 
phase  angle  as  defined  to  better  than  0.  02  or  0.  02%.  Thus,  the  quad¬ 
rature  of  the  phase-lock  amplifiers  is  checked  by  putting  bits  into  the 
signal  port  and  the  exclusive  OR  into  the  reference  channel.  The  90 
angles  as  read  out  from  the  phase-lock  amplifiers  were  found  to  be  good 
to  0.  05°.  The  fractional  quadrature  error  results  in  a  fraction  linewidth 
error  of  the  same  magnitude.  Therefore,  the  systematic  quadrature  error 
tends  to  cancel  in  the  subtractions  and  divisions  used  to  obtain  transport 
properties  from  linewidths. 

4.5  DATA 

Measurement  of  viscosity,  Prandtl  number  and  sound  speed  are 
presented  in  Tables  6-13.  The  gases  covered  are  argon,  nitrogen, 
hydrogen  and  air.  The  accuracy  of  the  argon  and  nitrogen  viscosity 
data  below  100  atm. is  0.  1%  in  viscosity  and  0.  5%  in  Prandtl  number. 

The  viscosity  data  in  this  range  is  within  0.  2%  of  the  oscillating  disk 
results. 43  The  remaining  data  have  not  been  compared  in  detail  with 
other  work  but  spot  checks  indicate  accuracies  better  than  1%. 

The  maximum  temperature-pressure  point  measured  was  405 
atm  at  600  K.  The  atmospheric  data  in  nitrogen  and  air  were  extended 
to  1000  K.  End  plate  resonances  prevented  data  taking  in  hydrogen  at 
elevated  temperatures. 


*wm  refers  to  the  frequency  of  the  fundamental  component  of  the  square 
wave. 


76 


Table  6.  High  temperature  measurements  in  nitrogen  (1  atm). 


Temperature 

°K 

Viscosity 

-4 

x  10 

Poise 

-1  -1 

g  sec  cm 

Prandtl 

Number 

Thermal 

Conductivity 

-4 

x  10 

Watts  cm  *  °K  * 

Specific 

Heat 

Ergs  g  1  °K  1 
(Ref.  33) 

Density  . 
-3 

g  cm 

(Ref.  26) 

299.0 

1.  789 

0.  714 

2.  619 

1.  040 

1. 1414 

382.  5 

2.  159 

0.  705 

3.  192 

1.  043 

0.  8922 

407.  1 

2.  283 

0.  708 

3.  368 

1.  045 

0.  8383 

420.  3 

2.  306 

0.  700 

3.  445 

1.  046 

0.  8120 

435.1 

2.  378 

0.  720 

3.  553 

1.  048 

0.  7843 

437.  2 

2.  374 

0.  697 

3.  568 

1.  048 

0.  7806 

440.  7 

2.  380 

0.  693 

3.  601 

1.  048 

0.  7744 

608.  9 

2.  984 

0.  695 

4.  615 

1.  076 

0.  5605 

762.  0 

3.  473 

0.  694 

5.  562 

1.  112 

0.  4479 

856.  6 

3.  751 

0.  699 

6.  096 

1.  136 

0. 3984 

950.  6 

4.  031 

0.  705 

6.  619 

1.  157 

0.  3590 

992.  3 

4.  119 

0.  697 

6.  893 

1.  166 

0.  3439 
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Table  7.  High  temperature  measurements  in  air  (1  atm). 


Viscosity 

Thermal 

Specific 

-4 

x  10 

C  onductivity 

Heat 

Density 

-4 

-1  o  -1 

-3 

Temperature 

Poise 

Prandtl 

x  10 

Ergs  g  K 

g  cm 

°K 

-1  -1 

g  sec  cm 

Number 

Watts  cm"1  °K-1 

(Ref.  33) 

(Ref.  27) 

350.  6 

2.  109 

0.  713 

2.  98 

1.  009 

398.  1 

2.  315 

0.  710 

3.  30 

1.  014 

404.  4 

2.  342 

0.  701 

3.  34 

1.  014 

0.  873 

425.  1 

2.  429 

0.  709 

3.47 

1.  017 

0.  830 

591.  8 

3.  064 

0.  709 

4.  53 

1.  049 

0.  596 

599.  2 

3.  104 

0.  709 

4.  50 

1.  051 

0.  589 

716.4 

3.  486 

0.  711 

5.29 

1.  079 

0.  493 

750.  3 

3.  593 

0.  711 

5.  49 

1.  087 

0.  470 

799.4 

3.  758 

0.  707 

5.  75 

1.  097 

0.  441 

805.  5 

3.  777 

0.  703 

5.  90 

1.  100 

0.  438 

895.  2 

4.  051 

0.  710 

6.  38 

1.  120 

0.  394 

954.  9 

4.  224 

0.  715 

6.  59 

1.  132 

0.  370 

380.  2 

2.  241 

0.  711 

3.  19 

1.  012 

0.  928 

401.  6 

2.  330 

0.  711 

3.  33 

1.  014 

0.  879 

448.  3 

2.  526 

0.  709 

3.  63 

1.  021 

0.  787 

519.  2 

2.  804 

0.  708 

4.  09 

1.  033 

0.  680 

609.  1 

3.  137 

0.  705 

4.  69 

1.  054 

0.  579 

650.  7 

3.  260 

0.  710 

4.  87 

1.  063 

0.  542 

699.  0 

3.  446 

0.710 

5.  20 

1.  074 

0.  505 

711.  6 

3.  467 

0.  710 

5.  25 

1.  077 

0.  496 

790.  1 

3.  726 

0.  703 

5.  69 

1.  097 

0.  447 

850.  6 

3.  909 

0.  703 

6.  07 

1.  110 

0.  415 

880.  2 

3.  987 

0.  697 

6.  39 

1.  117 

0.401 

956.  4 

4.  241  | 

0.  703 

6.  70 

1.  133 

0.  369 
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Table  8.  High  pressure  measurements  in  nitrogen. 


Viscosity 

Specific 

Thermal 

x  10" 

4 

Heat  Cp 

Specific 

Sound 

C  onductivity 

Temperature 

Pressure 

Poise 

P  randtl 

x  107 

Heat 

Speed 

x  10“4 

°K 

atm 

-1 

g  sec 

-1 

cm 

Number 

„  -1  o-l 

Ergs  g  K 

Ratio 

cm(p.sec)  1 

Watts  cm  1  °K  * 

294.  2 

1.  000 

176. 

4 

0.  713 

1.  040 

1.400 

0. 03498 

257.  3 

291.  7 

4.  940 

175. 

1 

0.  710 

1.  061 

1.402 

0. 03491 

261.  7 

289.  2 

5.  300 

174. 

7 

0.  712 

1.  062 

1.401 

0. 03478 

260.  6 

292.  2 

6.  700 

175. 

8 

0.  713 

1.  065 

1.408 

0. 03497 

250.  1 

295.  4 

57.  82 

187. 

4 

0.  724 

1.  128 

1.499 

0. 03619 

291.9 

295.  3 

66.  53 

189. 

6 

0.  728 

1.  142 

1.  513 

0. 03640 

297.4 

294.  8 

72.  86 

191. 

0 

0.  732 

1.  152 

1.523 

0. 03669 

300.6 

298.  9 

87.  62 

196. 

9 

0.  730 

1.  1 6 1 

1.  550 

0. 03723 

313.  2 

298.  8 

101.  2 

200. 

4 

0.  728 

1.  177 

1.  571 

0.  03767 

324.0 

298.  9 

102.  8 

201. 

1 

0.  728 

1.  779 

1.  573 

0. 03768 

325.  7 

298.  8 

108.  5 

202. 

1 

0.  733 

1.  190 

1.  582 

0. 03797 

328.  1 

299.  0 

112.  6 

204. 

0 

0.  729 

1.  196 

1.  589 

0. 03815 

334.  7 

298.9 

115.  4 

204. 

2 

0.  732 

1.  197 

1.  593 

0. 03824 

333.9 

299.0 

127.  9 

208. 

8 

0.  736 

1.  214 

1.  613 

0. 03874 

344.4 

298.  6 

137.  7 

209. 

9 

0.  739 

1.  227 

1.  624 

0. 03919 

348.  5 

299.  6 

159.  9 

219. 

5 

0.  742 

1.  240 

1.  674 

0.  04701 

366.  8 

297.  6 

240.  3 

244. 

0 

0.  747 

1.  309 

1.  831 

0. 04352 

426.4 
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Table  9.  High  pressure  measurements  in  air. 


Temperature 

°K 

Pressure 

atm 

Viscosity 

-4 

x  10 

Poise 

-1  -1 
g  sec  cm 

Prandtl 

Number 

Specific 

Heat  Cp 
x  107 

Ergs  g"1  °K_1 

Specific 

Heat 

Ratio 

Sound 

Speed 

/  i*l 

cm(|i  sec) 

Thermal 

Conductivity 
x  10'4 

Watts  cm  *  °K  * 

298.  6 

10.  01 

1.  886 

0.  711 

1.  012 

1.  417 

0. 03529 

2.  686 

299.  7 

21.42 

1.  916 

0.  710 

1.  020 

1.465 

0. 03550 

2.  749 

298.  8 

36.  30 

1.  937 

0.  702 

1.  034 

1.486 

0.  03568 

2.  855 

298.  1 

52.  7 

1.  968 

0.  722 

1.  087 

1.498 

0.  03587 

2.  960 

297.  8 

91.  5 

2.  051 

0.  739 

1.  148 

1.  554 

0. 03669 

3.  186 

298.  0 

101.  2 

2.  104 

0.  747 

1.  1 6 1 

1.  572 

0. 03699 

3.  269 

296.  8 

108.  2 

2.  119 

0.  752 

1.  171 

1.  592 

0. 03713 

3.  295 

297.  1 

116.  4 

2.  149 

0.  762 

1.  185 

1.  598 

0. 03740 

3.  339 

297.  2 

231.  1 

2.  175 

0.  762 

1.  194 

1.  638 

0,  03763 

3.  410 

297.4 

130.  5 

2.  198 

0.  768 

1.  202 

1.  609 

0. 03790 

3.  442 

298.  3 

152.  1 

2.995 

0.  787 

1.  232 

1.649 

0. 03881 

3.  592 

298.  4 

211.  8 

2.  603 

0.  840 

1.  302 

1.  704 

0. 04167 

3.  994 

298.  3 

299.  5 

2.  904 

0.  835 

1.  363 

1.  913 

0. 04665 

4.  743 

Table  10.  High  pressure  measurements  in  hydrogen. 


Viscosity 

Specific 

SBHfipjl 

-4 

x  10 

Heat  Cp 

Specific 

Sound 

Thermal 

Temperature 

Pressure 

Poise 

Prandtl 

x  107 

Heat 

Speed 

Conductivity 

atm 

-1  -1 
g  sec  cm 

Number 

Ergs  g”1  °K'1 

Ratio 

cm(M’Sec) 

Watts  cm”1  °K-1 

295.  15 

1. 441 

0.  884 

0.  703 

14.  311 

1.  405 

1301.  9 

17.  9 

295.  29 

3.  945 

0.  884 

0.  703 

14.  319 

1. 405 

1302.  8 

18.  0 

295. 26 

5.  556 

0.  884 

0.  704 

14.  328 

1.403 

1307.  8 

17.  9 

295.  25 

6.  610 

0.  884 

0.  702 

14.  328 

1.405 

1307.  9 

18.  0 

297.  16 

9.  120 

0.  883 

0.  702 

14.  336 

1.406 

1307.  4 

18.  0 

291.  08 

15.  65 

0.  882 

0.  702 

14.  315 

1.  406 

1314.  6 

17.  9 

291.  79 

22.  31 

0.  883 

0.  703 

14.  369 

1.  406 

1321.  8 

18.0 

296.  73 

24.66 

0.  884 

0.  704 

,  14.  377 

1.  407 

1331.  6 

18.  0 

297.  78 

33.23 

0.  883 

0.  697 

14.  398 

1.  407 

1342.  4 

18.2 

297.  57 

39.  80 

0.  885 

0.  702 

14.  414 

1.  407 

1326.  3 

18.2 

297.  47 

42.  17 

0.  884 

0.  702 

14.  418 

1.  408 

1327.  5 

18.2 

295. 26 

48.  30 

0.  885 

0.  702 

14.  439 

1.  409 

1350.  3 

18.  2 

297.  56 

56.  90 

0.  888 

0.  701 

14.  460 

1.  409 

1359.  6 

18.  3 

269-  31 

57.  80 

0.  890 

0.  702 

14.  460 

1. 410 

1366.  6 

18.  3 

296.  63 

63.  30 

0.  891 

0.  702 

14.  472 

1.410 

1365.  9 

18.4 

294.  60 

69.  00 

0.  894 

0.  703 

14.  484 

1. 411 

1368.  7 

18.4 

299.  19 

88.  80 

0.  908 

0.  702 

14.  526 

1.414 

1393.  7 

18.  8 

299. 15 

108.  8 

0.  912 

0.  703 

14.  550 

1.  420 

1411.  4 

18.  9 
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Table  11.  High  pressure  measurements  in  argon 


Temperature 

°K 

Pressure 

atm 

Viscosity 
x  10-4 

Poise 

-1  -1 
g  sec  cm 

Prandtl 

Number 

Specific 

Heat  Cp 

7 

x  10 

Ergs  g"1  °K"1 

Specific 

Heat 

Ratio 

Sound 

Speed 
cm(psec)  1 

Thermal 

C  onductivity 
Watts  cm  *  °K  * 

297.  2 

0.  500 

2.  256 

0.  689 

0. 5209 

1.  670 

. 003200 

1.  73 

297.  2 

1.  00 

2.  255 

0.  89 

0. 5216 

1.  670 

0.  03203 

1.  70 

299.  2 

20.  41 

2.  302 

0.  691 

0.  5459 

1.  728 

0.  03235 

no 

00 

H 

294.  8 

26.  80 

2.  295 

0. 690 

0. 5543 

1.749 

0.  03213 

1.  84 

296.  9 

28.  10 

2.  319 

0.  692 

0. 5668 

1.  746 

0.  03225 

1.  90 

299.  2 

40.  82 

2.  363 

0.  7,00 

0. 5703 

1.  791 

0. 03255 

1.  93 

295.  7 

46.  50 

2.  358 

0.  702 

0. 5734 

1.  811 

0.  03250 

1.  93 

294.  4 

67.  60 

2.  437 

0.  725 

0. 6119 

1.  873 

0.  03278 

3.  06 

294.  1 

81.  36 

2.  480 

0.  739 

0. 6408 

1.  947 

0. 03286 

2.  15 

296.  0 

100.  0 

2.  571 

0.  755 

0. 6589 

1.  971 

0. 03328 

2.  24 

295.  4 

103.  6 

2.  582 

0.  754 

0.  6621 

1.  986 

0. 03331 

2.  27 

294.  2 

110.  3 

2.  610 

0.  762 

0.  6712 

2.  024 

0. 03344 

2.  30 

294.  2 

117.4 

2.  645 

0.  776 

0. 6785 

2.  048 

0.  03361 

2.  31 

295.  2 

125.  1 

2.  700 

0.  777 

0. 6859 

2.  063 

0. 03389 

2.  38 

295.  0 

132.  3 

2.  724 

0.  776 

0. 6938 

2.  094 

0.  03430 

2.  43 

AEDC-TR-74-33 


Table  12.  High  pressure-temperature  measurements  in  argon. 


Pressure 

atm 

Viscosity 
x  10"4 

Poise 

-1  -1 
g  sec  cm 

Prandtl 

Number 

Specific 

Heat  Cp 
x  107 

Ergs  g’1  °K"1 

Specific 

Heat 

Ratio 

Sound 

Speed 
cm(ixsec)  * 

Thermal 

C  onductivity 
Watts  cm  *  °K  * 

425.  7 

20.42 

3.  075 

0.  732 

0. 5761 

1.  7061 

0. 03801 

2.42 

426.  1 

59.  70 

3.  141 

0.  725 

0. 5865 

1.  721 

0. 03942 

2.  54 

423.  1 

108.  2 

3.  565 

0.  778 

0. 5824 

1.  787 

0. 03997 

2.  67 

424.  8 

153.  2 

3.  350 

0.  737 

0. 6033 

1.  803 

0. 04273 

2.  74 

522.  5 

214.  8 

3.  511 

0.  737 

0. 6033 

1.  956 

0. 04372 

3.  01 

518.  7 

254.  6 

3.  984 

0.  786 

0. 6316 

1.  899 

0. 04491 

3.  37 

521.  5 

351.  1 

4.  119 

0.  744 

0. 6656 

1.  982 

0. 04982 

3.  84 

522.  6 

365.  8 

4.  211 

0.  762 

0. 6891 

2.  014 

0. 04995 

3.  81 

510.5 

415.  6 

4.  354 

0.  7.52 

0. 6895 

2.  263 

0.  04521 

3.  88 

510.  4 

25.  33 

3.  978 

0.  744 

0. 6702 

1.  685 

0. 04141 

2.  77 

507.  7 

154.  7 

3.  736 

0.  706 

0. 5182 

1.  829 

0.  04283 

3.  08 

507.  2 

221.  7 

3.  843 

0.  699 

0. 5823 

1.  932 

0.  04352 

3.  25 

511.3 

330.  1 

3.  596 

0.  612 

0. 5914 

2.  019 

0. 04578 

3.  52 

601.  9 

480.  1 

4.405 

0.  701 

0. 5998 

2.  196 

0. 04792 

3.  78 

602.  3 

220.  5 

4.  231 

0.  704 

0.  6021 

1.  834 

0. 05814 

3.  52 

600.  0 

318.  7 

4.  393 

0.  700 

0. 5863 

1.  895 

0.  06013 

3.  70 
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Table  13.  High  pressure-temperature  measurements  in  nitrogen. 


Viscosity 

Specific 

x  10'4 

• 

Heat  Cp 

Specific 

Sound 

Thermal 

Temperature 

Pressure 

Poise 

Prandtl 

x  107 

Heat 

Speed 

Conductivity 

°K 

atm 

-1  -1 
g  s  ec  cm 

Number 

Ergs  g"1  °K-1 

Ratio 

cm(psec)  1 

Watts  cm  ^  °K  * 

392.  6 

10.  29 

2.  207 

0.  714 

1.  050 

1.  412 

0.  04063 

3.  25 

421.  9 

51.  24 

2.  373 

0.  720 

1.  093 

1.  432 

0.  04312 

4.  14 

401.  2 

99.  94 

2.  370 

0.  725 

1.  115 

1.  474 

0.  04355 

3.  64 

422.  6 

152.  3 

2.  563 

0.  718 

1.  071 

1.  548 

0.  04727 

3.  51 

398.4 

214.  3 

2.  595 

0.  775 

1.  013 

1.  532 

0.  04763 

3.  39 

391.2 

238.  2 

2.  641 

0.  772 

1.  197 

1.  552 

0.  04843 

4.  09 

392.  6 

422.  9 

3.  144 

0.  781 

1.  175 

1.  653 

0.  05856 

4.  73 

501.  2 

11.  13 

2.  620 

0.  712 

1. 056 

1.  398 

0.  04576 

3.  89 

501.  6 

172.  2 

2.  817- 

0.  765 

1.  095 

1.  483 

0.  05104 

4.  02 

519.  4 

219.  7 

3.  057 

0.  790 

1.  167 

1.  896 

0.  06029 

4.  52 

511.  1 

521.  1 

3.  641 

0.  785 

1.  175 

1.  901 

0.  0658 

5.  45 

620.  5 

174.  9 

3.  177 

0.  725 

1.  061 

1.462 

0.  05617 

4.  65 

622.  37 

192.  9 

3.  224 

0.  782 

1.  070 

1.  797 

0.  06269 

4.41 

624.  0 

413.  5 

3.  498 

0.  775 

1.  162 

1.  671 

0.  06802 

5.24 

624.  1 

502.  3 

3.  139 

0.  773 

1.  165 

1.  696 

0.  07181 

4.  73 
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SECTION  V 
SUMMARY 

Acoustic  resonators  in  gases  provides  a  means  for  extremely 
precise  (1-0.  01%)  thermal  conductivity  and  viscosity  determinations. 

In  addition,  sound  speed  is  obtained  at  each  pressure,  volume,  tem¬ 
perature  point  so  that  specific  heat  (sound  speed)  and  relaxation  time 
data  are  also  obtained.  Therefore  a  measurement  with  the  acoustic 
resonator  technique  provides  more  information  than  separate  mea¬ 
surements  on  a  conventional  thermal  conductivity,  viscosity  plus  a 
Joule-Thomps on  apparatus.  Thus,  one  measurement  on  an  acoustic 
resonator  provides  as  much  data  as  three  conventional  measurements. 

In  addition,  the  data  is  self-consistent  and  potentially  one  to  two  orders 
of  magnitude  more  accurate  than  existing  techniques. 

Theoretical  expression  for  the  linewidth  and  center  frequency  of 
the  resonator  is  worked  out  for  a  nonideal  gas  with  relaxation  effects. 

Slip,  accommodation  and  acoustic  propagation  into  the  walls  are  in¬ 
cluded.  Slip  corrections  to  viscosity  and  accommodation  corrections 
for  thermal  conductivity  are  ~0.  5%  at  300°K  (1  atm)  and  increase  to 
~0.  8%  at  1000  K  (1  atm)  for  a  1  cm^  volume  cavity.  However,  these 
corrections  become  smaller  directly  with  the  square  root  of  density  and 
inversely  with  the  square  root  of  a  typical  cavity  length.  The  losses 
due  to  propagation  into  the  cavity  walls  increase  from  0. 1%  in  linewidth 
at  ~1  amagat  to  ~20%  in  linewidth  at  ~40Q  amagats.  Methods  for  calcu¬ 
lating  these  quantities  from  measured  quantities  are  detailed  in  the  text. 

Transducers  for  introducing  sound  into  a  resonator  and  sensing 
the  resulting  fields  also  contribute  to  the  linewidth  and  center  frequency. 
Methods  for  calculating  the  contribution  of  the  transducers  to  the  line- 
width  are  given.  The  data  reduction  techniques  for  obtaining  the  linewidth 
from  measurements  with  superimposed  electrical  feedthrough  and  mode 
coupling  signals  are  discussed.  These  transducer  associated  effects 
amounted  to  as  much  as  50%  of  the  measured  linewidths  at  high  pressures. 
However,  with  a  combination  of  calculation  and  self-consistent  iteration 
the  transducer  effects  can  be  handled  to  a  precision  consistent  with  ~0.  1% 
in  viscosity. 

Mechanical  details  of  the  transducers  and  cavities  used  in  the 
present  work  are  given.  Mica  dielectric  electrostatic  transducers 
coupled  to  the  cavity  with  coupling  holes  (~0.  2  cm  long,  0.  005  cm  radius) 
were  developed  for  use  up  to  1000°K.  The  tradeoffs  between  coupling 
corrections  and  allowable  signal  levels  as  a  function  of  coupling  hole 
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size  are  elaborated.  It  is  concluded  that  an  order  of  magnitude  re¬ 
duction  of  coupling  hole  correction  beyond  that  obtained  in  the  present 
work  is  feasible  by  specializing  the  hole  size  for  narrow  pvt  ranges. 
Difficulties  were  encountered  in  running  the  mica  dielectric  trans¬ 
ducers  at  temperatures  above  600  K  for  times  greater  than  ~8  hr. 
Therefore,  electrostatic  transducers  with  thick  (0.002")  quartz 
dielectric  plates  or  conduction  of  the  sound  from  regions  below  600  K 
with  acoustic  wave  guides  are  suggested  for  future  high  temperature 
work.  It  should  be  noted  that  the  acoustic  wave  guide  option  could 
extend  cavity  measurements  to  about  2000  K. 

Steel  and  invar  cavities  with  volumes  ranging  between  1  and  4 
cm^  used  in  the  present  work  are  described  in  the  text.  The  design 
tradeoff  of  oxidation  resistance  and  acoustic  impedance  versus  ease 
of  fabrication  and  procurement  is  discussed.  The  greatest  problem 
with  the  cavity  design  is  flexural  vibrations  in  the  end  plate.  Thus, 
high  modulus  materials  are  optimum  for  end  plate  designs.  It  is  also 
pointed  out  that  cavities  with  volumes  of  up  to  10  cm^  are  optimum 
for  work  at  high  temperatures  or  in  hydrogen  and  helium. 

Electronic  systems  for  the  accurate  measurement  of  phase  and 
amplitude  versus  frequency  are  described.  A  frequency  synthesizer 
and  a  quadrature  phase-lock  amplifier  are  described  which  allow  0.  01% 
linewidth  determinations  to  be  made.  A  digital  phase  shift  synthesizer 
is  described  along  with  a  modulation  system  capable  of  0.  001%  line- 
width  determinations.  However,  the  modulation  system  still  requires 
a  line  sweep  determination  to  eliminate  feedthrough  and  correct  for 
Lorentzian  line  shapes.  The  electronics  to  drive  and  bias  the  trans¬ 
ducers  is  also  described. 

Finally,  data  taken  with  the  systems  described  above  are  pre¬ 
sented.  Accuracies  of  0.  1%  in  viscosity  and  0.  5%  in  Prandtl  number 
are  given  for  hydrogen,  nitrogen,  air  and  argon  at  300  K  at  pressures 
up  to  about  150  atm.  Data  in  argon  and  nitrogen  are  presented  at  pres¬ 
sures  up  to  460  atm  and  temperatures  up  to  600°K.  In  addition,  mea¬ 
surements  in  nitrogen  and  air  at  temperatures  up  to  1000  K  and  at 
atmospheric  pressure  are  given.  These  last  two  sets  of  data  have  ac¬ 
curacies  approaching  0.  5%  in  viscosity  and  1%  in  thermal  conductivity 
on  the  high  temperature  and  high  pressure  ends. 
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APPENDIX  A 

TRAVELING  WAVES  IN  A  CYLINDRICAL  WAVE  GUIDE 


PHYSICAL  VARIABLES 

The  cavity  is  a  cylindrical  wave  guide  with  end  wall  mounted 
perpendicular  to  its  axis.  The  linewidth  and  center  frequency  expres¬ 
sions  can  be  obtained  by  first  making  wave  functions  which  satisfy  the 
wave  equations  and  the  boundary  equations  on  the  side  walls.  These 
waves  are  then  allowed  to  reflect  back  and  forth  between  the  end  walls 
of  the  cavity.  The  sum  of  all  these  reflections  is  the  desired  field  in 
the  cavity. 

The  wave  equations  are  expressed  in  terms  of  the  potentials  Hj, 
H2  and  H3  (see  Table  A.  1).  However,  the  boundary  conditions  apply  to 
the  pressure,  temperature  and  flow  velocity.  The  flow  velocity  is  re¬ 
lated  to  the  potentials  in  the  usual  way* »  ^ 

v  =  V  (HL  +  H2)  +  VxH3.  (A.  1) 


The  second  term  in  (A.  1)  can  be  simplified  if  only  axisymmetric  modes 
are  considered.  The  contribution  to  the  particle  velocity  due  to  the 
viscous  wave  is 


r0 

A 

z 

« 

9 

9 

9 

9  r 

90 

9  z 

H 

r 

He 

H 

z 

(A.  2) 


The  caret  indicates  unit  vectors  in  the  radial  (r),  azimuthal  r^  ,  and  axial 
(z)  direction.  The  axisymmetric  condition  means  that  the  azimuthal  com¬ 
ponent  of  the  particle  velocity  must  vanish.  Thus,  Hr  =  Hz  =  0.  In  addi¬ 
tion,  for  simplicity  we  modify  the  nomenclature  so  that  is  a  scalar 
satisfying 


9 


9  r  9  z  ’ 


(A.  3) 


Notice  that  this  new  potential  (H3)  satisfies  the  viscous  wave  equation.  In 
addition,  no  spurious  solutions  are  generated  by  the  transformation  in 


93 


AEDC-TR-74-33 


Eq.  (A.  3)  since  only  the  differentiated  form  appears  in  the  particle  ve¬ 
locity.  The  particle  velocity  in  the  radial  (r )  and  axial  (z  )directions 
are  from  (A.  1)  and  A.  3) 


r  = 
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8  r 


H1  +  Hz  + 


82H, 


8  z 
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(A.  4) 


(A.  5) 


The  acoustic  temperature  amplitude  T  is  given  in  terms  of  the 
quantity 


6  = 


— T  =  a  H  +  a  H_ 
7-1  11  2  2 


(A.  6) 


where  (3  is  the  thermal  expression  coefficient,  and  y  the  equilibrium 
specific  heat  ratio,  aj  and  are  given  in  Table  A.  1.  Equation  (A.  6) 
is  identical  to  (A.  23)  of  Ref.  3.  However,  a^  and  a ^  can  be  simplified  to 


ai  1  C 
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(A.  7) 


v 


v 


to  better  than  1  pt  in  10  for  the  conditions  of  the  present  work.  Terms 
of  the  order  T  have  been  neglected  compared  to  1  in  and  terms  of  the 


order  of  ^  ~  J  have  been  neglected  compared  to  1  in  a^.  The  approxi- 
2  2 


mation 


-ft) 
“ft) 
ft) 
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<<  1  is  employed  many  times  below.  Therefore,  values 


k2  \2 

of  I-:; — I  are  given  in  Table  A.  2  for  several  pressures  and  temperatures 


in  nitrogen.  I  ~r~~  I  gets  smaller  directly  as  density  over  thermal  con- 

\kz/ 

ductivity  are  increased  and  as  frequency  is  decreased.  Experimentally, 
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Table  A.  2. 


(Ratio  thermal  to  acoustic  wavelength)  for  nitrogen 

* 

at  several  pressures  and  temperatures. 


Temperature 


o 


K 


300 

1000 

1.  01  x  10"3 

3.  01  x  10'"3 

1.  18  x 10-4 

4.  07  x  10“4 

0.  71  x  10“4 

2.  1  x  10"4 
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this  fact  makes  high  pressures  and  low  frequencies  favorable. 


Finally,  an  expression  for  acoustic  pressure  amplitude  (p)  in 
terms  of  potentials  is  obtained  from  the  equations  of  state  and  con¬ 
servation  of  energy.  The  conservation  of  energy  is  (a  corrected 
version  of)  Eq.  (A.  16)  of  Ref.  1 


V  0 


+ 


(A.  8) 


Differentiating  the  equation  of  state  with  respect  to  time  as  using  (A.  8) 


to  eliminate 


the  equation  of  state  becomes 
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Here  c  is  the  equilibrium  sound  speed  and  pQ  the  ambient  density.  The 
energy  equation,  together  with  the  simple  harmonic  dependence  of  the 

— i  u>t 

pressure  amplitude  (e  ),  allow  (A.  9)  to  be  transformed  into- 
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[v2  fl  +  7  C  fl]  . 


(A.  10) 


This  can  be  further  simplified  using  (A.  6)  to  eliminate  9 
expression  for  pressure  becomes 


The  final 


p  =  p  c  a 


*  cT  7Hi  °fe)  ' 


(A.  11) 


*More  mundane  considerations  of  spurious  resonances  in  the  cavity  walls 
favor  small  cavities.  Also,  at  very  high  densities  acoustic  propagation 
into  the  walls  becomes  important. 

Note  we  have  used  C  -Cv  =  Cp-Cv  since  the  contribution  of  internal 
modes  (vibration,  etc.  j  is  the  same  for  Cp  and  Cv  and  is  independent 
of  pressure  for  gases. 
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APPLICATION  OF  BOUNDARY  CONDITIONS 


The  wave  equations  in  Table  A.  1  may  be  separated  in  cylindrical 
coordinates.  In  fact,  for  the  cylindrical  wave  guide  a  complete  set  of 
orthogonal  solutions  can  be  found  exactly.  This  set  of  wave  functions 
are  then  used  to  obtain  an  extremely  precise  solution  to  the  cavity  prob¬ 
lem.  In  addition,  this  complete  orthogonal  set  of  cylindrical  wave  guide 
solutions  serves  as  the  basis  of  transducer  coupling  corrections.  There¬ 
fore,  we  have  gone  into  the  cylindrical  wave  guide  problem  in  great  detail. 


The  wave  equations  in  Table  A.  1  can  be  separated  in  the  usual  way 
into  the  solutions  of  the  type 


(A.  12) 


The  z  dependent  terms  in  the  viscous  and  thermal  waves  can  be  suppressed 
for  the  cylindrical  wave  guide  problem  since  the  wave  numbers  involved 
have  large  imaginary  parts.  Thus,  they  are  rapidly  damped  any  distance 
greater  than  ~0.  1  mm  from  the  wall.  The  conditions  of  separation  require 


that 


(A.  13) 


be  satisfied  by  the  separation  constants 


and 


7 


(3) 

P 
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The  three  equations  (A.  13)  can  be  used  to  eliminate  all  the  constants 

(l)2 

of  separation  in  terms  of  k  and  the  known  wave  numbers  k  ,  k_ 

and  k3>  1  1  L 


Thus, 


It  is  already  known  that  the  losses  have  a 

v  <l).  V1’ 

\  P  P- 

is  of  the  order  of  the  cavity  radius  for  the  first  mode 


20%)  effect  on  the  wave  numbers 

*r,l) 


:  cl 

)°f 


small"  (less  than 
the  acoustic  mode. 


(p  =  0)  in  Eq.  (A.  13).  The  higher  p'th  order  modes  are  radial  modes 
which  travel  even  slower  than  the  zero  order  mode.  That  is,  for  p>  0 


k(1» 

P 


2 


Therefore,  the  approximations 


7 


P 


(A.  14) 


(A.  15) 


(A.  16) 


can  be  made.  These  approximations  have  been  explored  above  (see 
Table  A.  2).  However,  the  theory  can  be  worked  through  with  the  ap¬ 
proximate  radial  wave  numbers  and  a  correction  made  using  these 

answers  which  will  be  good  to  ~ 

The  boundary  conditions  are  determined  by  the  physical  processes 
which  are  taking  place.  Energy  and  momentum  exchange  at  the  wall  must 
be  taken  into  account  at  low  densities  (~1  amagat)  and  high  precisions 
(~1  pt  in  10"4  in  linewidth).  Acoustic  transmission  into  the  walls  of  the 
cavity  becomes  important  at  elevated  densities  (~100  amagat).  These 
phenomena  cause  the  fluctuating  components  of  the  particle  velocity  and 
temperature  to  be  different  from  zero  on  the  cavity  walls. 
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The  radial  component  of  the  particle  velocity  at  the  wall  (rw) 
differs  from  zero  mainly  due  to  the  acoustic  propagation  into  the- 
walls  of  the  wave  guide.*  This  problem  is  quite  complex  to  solve 
in  general. 4  However,  the  dependence  of  these  losses  and  pc  can 
be  shown.  The  frequency  or  p  c  dependence  of  the  acoustic  transmis¬ 
sion  losses  can  be  used  to  determine  the  corresponding  increment  in 
linewidth  experimentally.  The  transmission  losses  can  be  estimated 
to  within  an  order  of  magnitude^  by  assigning  a  modulus  Y  =  pw  cw^ 
to  the  walls.  The  modulus  is  defined  by 


r  =  2  it  iu  a 
1  w 


p  c 
w  w 


(A.  17) 


Here,  a  is  the  cavity  radius,  pw  the  wall  density  and  cw  the  sound  speed 
based  on  the  modulus  Y.  The  pressure  can  be  eliminated  using  (A.  11), 
giving  us 


w 


H 


1  P 


-P-9 

c 

w  w 


(A.  18) 


as  the  appropriate  boundary  condition. 


The  possibility  of  waves  reentering  the  cavity  once  they  have  been 
transmitted  into  the  wall  has  been  neglected.  The  contribution  of  such 


waves  is  of  the  order  of 


The  effective  modulus 


of  the  wall  p  w  cw  will  depend  on  the  type  of  deformation  induced  on  the  wall 
by  the  pressure  in  the  gas.  However,  on  thick-walled  cavities  this  will 
always  lie  somewhere  between  the  longitudinal  and  the  flexural  modulus 
for  the  wall  material.  The  axial  velocity  at  the  wall  is  not  affected  much 
by  elastic  motion  of  the  wall.  This  is  because  the  sound  speed  is  much 
faster  in  the  solid  than  in  the  gas.  Thus,  the  pressure  in  the  wall  tends 
to  be  uniform  over  the  length  of  the  cavity. 


However,  the  axial  velocity  is  different  from  zero  at  low  pressures 
due  to  slip.  The  axial  velocity  at  the  wall  is  given  by^ 


z 


i 

V 


8  z 

9  r  * 


(A.  19) 


*The  end  wall  losses  due  to  acoustic  transmission  affect  z  mainly. 
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The  accommodation  length  i  v  is  a  few  mean  free  paths  in  the  gas. 
The  particle  velocity  due  to  slip  is  very  small  and  is  only  important 
at  densities  ~1  amagat  and  below. 


Similarly,  the  temperature  at  the  wall  depends  on  thermal 
accommodation  between  the  gas  and  the  wall.  In  addition,  thermal 
waves  propagate  into  the  wall  and  cause  an  additional  wall  tempera¬ 
ture  variation.  Since  the  temperature  fluctuations  do  not  penetrate 
the  wall  very  far  the  one  dimensional  heat  equation^ 


82T 


8  r 


P  w  ^vw  8  T 

V  8t 

w 


(A.  20) 


describes  the  heat  conduction.  The  solution  to  (A.  20)  which  is  zero 

for  r — ►<»  and  T  e  at  the  walls  is 
w 


T 


esq} 


P  c 

W  vw 


2  X 


w 


(A.  21) 


Thus,  the  heat  flux  into  the  wall  Q  is 

w 


\ 

— ^  T  .  (A.  22) 

1  D  w 

The  heat  flux  Qw  to  the  wall  through  the  thin  layer  of  gas  with  a  few  mean 
free  paths  (1  t)  of  the  wall  is 


Q  = 

w 


w 


8  T 
8  r 


up  C  X. 
w  vw  w 


r  =  a 


T  = 

w 


Q 

w 


(T  -  T  )  = 
w 


(A.  23) 


where  T  is  the  temperature  on  the  outer  edge  of  the  accommodation  layer 
and  f  t  is  the  accommodation  length.  Equation  (A.  23)  is  the  boundary  con¬ 
dition  which  must  be  satisfied  by  the  hydrodynamic  waves. 


The  wall  temperature  can  be  eliminated  between  (A.  22)  and  (A.  23) 


T 


(A.  24) 
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Table  A.  4.  Simultaneous  equations  for  derived  from  temperature 

and  particle  velocity  boundary  conditions. 
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The  boundary  condition  on  T  thus  goes  to  T  =  0  as  the  ability  of  the 

wall  to  absorb  heat  becomes  very  large,  that  is,  as  r— —  1  and  i 

K  D  t 

w 

go  to  zero.  The  accommodation  length  i  ^  is  a  few  mean  free  paths 
long  (~10"5  cm  at  S.  T.  P.  ).  However,  as  the  pressure  is  dropped  i  t 
rapidly  approaches  to  the  thermal  wavelength  and  can  no  longer  be 
neglected.  The  thermal  diffusion  length  becomes  significant  at  elevated 
densities  (p  ~  500  amagats).  However,  this  does  not  become  clear  until 
the  boundary  conditions  are  applied. 

The  potentials  from  our  trial  solution  (A.  12)  can  now  be  put 
into  the  expressions  for  particle  velocity  and  temperature  (A.  4-A.  6) 
to  obtain  an  equation  for  the  allowed  values  of  radial  wave  number 
The  other  constants  are  already  known  from  the  separation 
conditions  (A.  13  and  following  discussion).  The  particle  velocity  and 
temperature  at  the  wall  are  given  by  the  equation  in  Table  A.  3  for  each 
value  of  p.  Here  we  have  used  the  formulas  relating?  J'0  and  J'j  to  JQ 
and  J  to  eliminate  the  derivatives  of  the  Bessel  functions.  Applying  the 
boundary  conditions  (A.  18),  (A.  19)  and  (A.  24)  to  (A.  25)  we  get  a  set  of 
simultaneous  equations  for  the  Ap'n),  shown  in  Table  A.  4. 

Notice  that  the  complex  parts  of  (A.  18)  have  been  dropped  (the 

complex  radial  wave  number  7p(l)  has  been  replaced  with  a  real  radial 

wave  number  7  ) 

P 


d 

s 


2  TT 


«r 
_ o 

c 

8 


(A.  27) 


used  in  the  boundary  condition  for  i  (Eq.  A.  18).  The  condition  that 
(A.  26)  can  be  solved  for  the  Ap  leads  in  the  usual  way  to  the  dispersion 
relation  for  YpU),  which  is  given  in  Table  A.  5.  This  relation  is  very 
complex  as  it  stands.  The  approximations  of  (A.  15)  and  (A.  16)  go  a  long 
way  towards  making  it  possible  to  solve  the  dispersion  relation  analytically. 
In  addition,  the  Bessel  functions  in  Eq.  (A.  27)  can  be  eliminated  if  rela¬ 
tions  such  as 


J’ 


o 


(A.  29) 


can  be  found. 
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Table  A.  5.  Dispersion  relation  for  acoustic  waves  in  a  cylindrical  wave  guide 
with  slip  accommodation  and  acoustic  transmission  into  the  walls. 
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CYLINDRICAL  WAVE  GUIDE  DISPERSION  RELATION 

The  classical  approximations  for  the  Bessel  functions1  with  the 

arguments  |y^2^aj  and  |y^3^aj  are  based  on  the  large  imaginary 

(2)  (3)  2 

parts  of  y^  '  and  y^  .  Recalling  (A.  15)  and  from  Table  A.  1 


>(2)  s  ^  Uc2i~(1  +  i) io3 


(A.  30) 


where  ]  lc^  1  is  of  the  order  of  10^  at  about  1  atm  and  increases  with  pres¬ 
sure.  Thus,  the  asymptotic  expansions  of  Hankel^  should  give  good  ap¬ 
proximations  here.  Using  9.  2.  5  and  9.  2.  11  of  Ref.  7,  we  find  letting 

(2)  (3) 

z  =  y  a  or  y  a 

P  P 


J^z) 

JIT) 

o 


tan  (z 


(A.  31) 


and  expanding  the  tangent  in  exponentials 

tan  (z  -  |)  =  i  o(e'2  I  Zl).  (A.  32) 

Thus,  the  approximation  (A.  31)  is  good  to  parts  in  10^  and  (A.  32)  in  parts 
in  e-^000.  Thus,  the  desired  Pj  for  the  thermal  and  viscous  waves  in 
(A. 29) are 


(A.  33) 


The  approximation  for  the  acoustic  mode  |yp^a  j  follows  from 

the  fact  that  the  acoustic  mode  is  not  changed  much  by  losses.  Therefore 

we  try  an  expansion^  of  J  / y  ^a  ]  in  a  series  consisting  of  the  solutions 
/  (1 )  \  •  1  \  P  f 

ly  a  j  to  the  cylindrical  wave  guide  problem.  Thus,  from  standard 
formulas  for  expansions  in  Bessel  functions'^ 
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J  (t  (1  >r )  =  5^  A  J  (7  r) 

°\  p  /  n  o  n 


(A.  34) 


where 


An  -  T  2 

t  J<X*>  '° 


/*  a 

l  /  r  hnr)  ^  '  (A'35) 


The  integration  indicated  in  (A.  35)  may  be  carried  out.  In  addition, 
we  can  use  the  boundary  condition  for  lossless  waves  J1  (y^a)  =  0. 


The  expansion  (A.  34)  then  becomes,  after  solving  for  J' 
(which  cornea  from  the  integral  in  (A.  35)) 


The  term  in  brackets  is  always  less  than 


1 

2 


~  10  and  decreases 


with  increasing  pressure  and  frequency.  In  addition,  the  summation 
in  (A.  25)  rapidly  approaches  1  as  the  radial  mode  number  p  becomes 
larger  than  p  =  1. 
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EXPRESSIONS  FOR  7  (1  *  AND  k  (1  * 

P  P 


The  dispersion  relation  (A.  28)  can  now  be  simplified  to 


The  dispersion  relation  solved  for  p  j  in  a  way  which  makes  all  the 
corrections  due  to  nonideal  boundary  condition  clear,  becomes 


7 


P 


+ 


(A.  39) 


-ia 

s 


.PS.. 

c 

w  w 


viscous 


The  first  term  on  the  right-hand  side  of  (A.  39)  consists  of  the 

(1)  5 

"losses”.  7  i  ^  can  be  neglected  to  parts  in  10  here.  The 
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effect  of  the  momentum  accommodation  length  is  to  increase  the  viscosity  ' 
by  a  factor 


1 


(A.  40) 


This  factor  can  easily  be  1%  in  a  small  cavity  at  ~1  atm.  However,  it 
is  at  most  0.  1%  for  the  conditions  of  the  present  work. 


The  second  term  is  the  thermal  conduction  term  plus  a  second 
order  term  involving  viscous  and  thermal  effects.  The  accommodation 

(D. 


1  -  7. 


effects  in  this  second  order  term 


1-7. 


(2), 


can  be  dropped  to  parts 


5 

in  1 0  in  thermal  conductivity.  The  effect  of  thermal  capacity  of  the 
wall  and  thermal  accommodation  at  the  surface  of  the  wall  is  to  increase 
the  thermal  conductivity  by  the  factor 


1 


(A.  41) 


The  effect  of  thermal  accommodation  f  t  is  important  only  at  the  lower 
pressures  as  indicated  above  for  momentum  accommodation.  The 
thermal  length  f  ^  of  the  wall  is  small  (parts  in  10“^)  at  about  1  amagat. 
However,  at  500  amagats  it  rises  to  about  1  pt  in  10“^. 


The  final  term  in  (A.  39)  is  the  loss  to  the  side  wall  by  acoustic 
transmission.  This  relation  ignores  any  frequency  or  wavelength 
behavior  due  to  the  geometry  of  the  cavity  side  walls.  However,  it 
shows  the  major  frequency  and  p  c  dependence  of  linewidth.  This  can 
be  used  in  practice  to  compensate  the  linewidth  measurements  for 
acoustic  transmission  losses  at  elevated  densities  as  described  in  the 
text  of  this  report.** 


As  is  seen  below 


**Notice  that  y 

P  v 


the  same  way! 


andlj_  .  /T. 

a2  V  0 

so  that  accommodation  can  be  handled 
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Since  we  can  easily  restore  the  slip,  accommodation  and  thermal 
capacity  corrections  by  modifying  viscosity  with  (A.  40)  and  thermal  con¬ 
ductivity  with  (A.  41),  we  will  no  longer  carry  these  terms  in  our  algebra. 
The  dispe-rsion  relation  then  becomes 


(1)  2 
7-7  =  — 

P  pa 


2  P*  2P2  -  ■  Po  k 


(ir  r  3 

’  „  (3)' 


(1) 


+  k!  ~2^-1)  +  kl  <7_1)  ZJT2 


(3)  k. 


*"  (l)2  2 

1  *  2 

t5^n  (1)  2 

7  -  7_ 

n  p 


+  i  a  7  k 


8  '  1  p  C 

w  w 


(A.  42) 


In  fact,  many  of  the  second  order  terms  can  be  considered  corrections 
to  the  viscosity  and  thermal  conductivity  in  the  manner  discussed  above 
for  accommodation  and  slip.  The  total  correction  factors  for  viscosity 
and  thermal  conductivity  are  shown  in  Table  A.  6.  The  remaining  ex¬ 
pression  for  the  radial  wave  number  is 


7 


(1) 

P 


-7 


P 


2 

a 


,2  ,2 

-  kl  -  kl 

+  (7-1)  T -  +  (7-1)  -T — 

k2  k3 


+  iff 

s 


pc 


c 

w  w 


(A.  44) 


The  expression  (A.  44)  is  the  basis  of  the  first  order  wave  number  ex¬ 
pressions.  The  true  viscosities  and  thermal  conductivities  must  be 
corrected  by  the  factors  in  Table  A.  6. in  order  to  obtain  wave  numbers 
correct  to  second  order  and  for  nonideal  boundary  conditions.  Typically, 
at  300  K,  1  atm  and  f01  ~  50  kHz  the  correction  factors  amount  to  ~0.  5% 
in  viscosity  and  decreases  with  increasing  pressure. 


110 


Table  A.  6.  Second  order  correction  factors  for  nonideal  acoustic  wave  guide 
in  p'th  radial  mode. 


Multiply  True  by  to  get  corrected  wave  numbers 


Thermal 
C  onductivity 
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The  wave  numbers  for  traveling  acoustic  waves  in  a  cylindrical 
wave  guide  are  to  the  first  order 


(1)  2  .  1  -  i  .  2  .  .  .  pc 

=  y  +  -  k  e  +  l  a  yk  — c — 


?  "T  ~sp  s  1  P  c 

c  w  w 


(A.  45) 


and  from  Eq.  (A.  1 3) 


(1)  .2  2  1  -  i  ,  2  .  pc  /a  4 

=  kl  “  7p  ~““k’  1  yK  ~  ~~  (A.  46) 


2  1  sp  s  1  1  p  c 

c  w  w 


where 


ft)- 


M- 

cop 


+  (7  -  1) 


VWPCT 


(A.  47) 


The  imaginary  part  of  kp  is  the  absorption  coefficient  for  pth  mode 
acoustic  wave  propagating  down  a  cylindrical  wave  guide.  The  real 
part  of  kp^'1 )  is  the  wave  number  (frequency/propagation  speed)  for 
such  waves.  The  y  (1 )  form  an  infinite  set  of  allowed  complex  wave 
numbers  closely  following  the  allowed  real  wave  numbers  of  the  loss¬ 
less  wave  guide  problem.  These  real  wave  numbers  are  the  successive 
roots  of  the  first  order  Bessel  functions  and  satisfy  the  equation 


J1  (7pa)  =  °  •  (A#48) 

We  can  now  make  an  orthogonal  set  of  wave  functions  out  of  the  acoustic 
mode  wave  functions  from  which  all  acoustic  disturbances  in  a  cylindrical 
wave  guide  can  be  made. 


THE  EFFECTS  OF  THE  SOUND  SOURCE 

/  •  v 

The  amplitude  of  the  modes  A  1  is  determined  by  the  velocity 
versus  radial  distance  on  a  source  located  at  z  =  0  in  the  wave  guide. 


The  functions  <b 

P 


satisfy  the  boundary  conditions 


on  the  side  wall  already.  Thus,  if  we  make  a  sum  of  which  satisfy  the 
boundary  conditions  on  the  end  wall  at  the  source  we  get  solutions  which 
satisfy  the  wave  equations  and  boundary  conditions  everywhere.  Transducers 
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which  cause  z  to  vary  only  with  r  at  the  source  are  used  in  the  present 
work.  The  expansion  of  z  Q  in  0  can  be  accomplished  to  give 

■T* 


z  = 


G  6 
P  P 


(A.  49) 


Here,  we  have  not  assumed  that  the  0  are  orthogonal.  They  are  not 
because  7p(l)  has  imaginary  parts.  However,  in  practice,  linear  sums 
of  the  <J>p  can  be  found  which  are  orthogonal  and  these  functions  used  to 
uniquely  calculate  the  Gp  of  Eq.  (A.  49).  This  fine  point  is  not  important 
for  the  cavity  because  as  already  noted  the  are  almost  orthogonal. 
However,  for  the  coupling  hole  problem  (a  ~  «£  or  ^3)  orthogonalization 
is  necessary  as  is  discussed  in  Appendix  C. 


Transducers  at  z  =  0  whose  particle  velocity  matches  the  particle 
velocity  of  the  pth  mode  exactly  will  generate  only  0 p.  Thus,  from  a 
physical  point  of  view  we  can  insist  that  each  </)p  satisfy  the  boundary 
conditions  on  the  source  wall  separately.  Plane  viscous  and  thermal 
wave  propagation  in  the  z  direction  are  required  to  satisfy  the  boundary 
conditions  at  the  wall  in  the  presence  of  the  acoustic  wave  <p p.  Thus  our 
trial  solution  (A.  12)  becomes 


(A.  50) 


Note  that  we  have  not  used  the  radial  viscous  and  acoustic  modes  since 

J  (7  ^  r  )  vanish  over  most  of  the  end  wall  (for  r  >  ~  ).  More  is 

o  Vp  /  K 


P  / 
said  about  this  below. 


Analogous  to  Eqs.  (A.  18),  (A.  19),  (A.  20)  the  boundary  conditions 
on  the  source  wall  for  the  pth  radial  mode  are  set  forth  in  Table  A.  7. 
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The  differences  in  the  end  wall  and  side  wall  conditions  are  due  to 
interchange  of  r  and  z  as  the  tangential  and  normal  velocity  com¬ 
ponents.  In  addition,  the  normal  velocity  is  specified  by  transducer 
design  in  the  coupling  coefficient  Gp. 

The  velocity  components  in  the  gas  just  outside  the  accom¬ 
modation  lengths  are  obtained  by  substituting  our  trial  solution  (A.  12) 

(with  J  ^7  ^  r  )  ,  i  =  2,  3  suppressed)  into  Eqs.  (A.  4),  (A.  5)  and 

(A.  6)  and  evaluating  them  at  z  =  0.  The  boundary  conditions  in  the 
gas  become 


(A.  52) 


Finally,  setting  (A.  52)  and  (A.  51)  equal  to  one  another  a  set  of  simul¬ 
taneous  equations  for  A  ^ ,  B  ^  and  B  ^  are  obtained 

P  P  P 
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(1) 


The  amplitude  of  the  wave  excited  by  the  source  transducer 
can  be  found  by  solving  (A.  5  3)  by  Cramer's  rule.  After  neglecting 

ai  -6 

—  (~10  )  compared  to  an  expression  formally  equivalent  to  Fritsche's 

a2 

equation  (33)  (Ref.  A.  1) 


(1) 


-iG 


k  +  h 


(A.  54) 


is  obtained.  However,  h,  including  nonideal  boundary  conditions,  is 


h  = 


1  -  ik  (1)i  7 

_ P  v  _ E 


1  -  ik  <3>i 


ID2 


p 


W 


k<3>  ) 

p  p  D  ‘/ 


(A.  55) 


and  k  ^  contains  all  the  second  order  terms  discussed  in  the  last  para- 

P  (1 ) 
graph.  The  most  significant  point  about  (A.  54)  is  that  every  time  k  —  0 

(1 )  ^ 

(recal  =  (A.  46))  the  amplitude  A^  gets  very  large.  Thus,  frequencies 

where  =  l7p|  are  radial  resonances.  The  coupling  coefficient  is  not 

strongly  dependent  on  frequencies  because  we  had  to  design  a  "broadband" 
transducer  for  the  application  under  study.  However,  multiple  reflections 
between  opposing  walls  will  pull  the  resonance  predicted  by-  (A.  54).  In 

...  -  2 

addition,  the  radial  resonance  goes  as  k  ^  ^  in  a  cavity ‘which  is  easier 

(1)  P 

to  calculate  than  the  k  dependence  in  an  infinite  wave  guide.  There¬ 
fore,  the  multiple  refaction  problem  in  a  cavity  with  two  end  walls  will 
be  done  before  the  modes  are  identified  in  detail. 

MULTIPLE  REFLECTIONS  AND  CAVITY  RESONANCES 

Walls  at  the  end  of  the  wave  guide  cause  the  waves  to  reflect  back 
and  forth.  The  total  field  at  any  point  is  the  sum  of  all  the  multiple  re¬ 
flections.  The  basic  traveling  acoustic  mode  wave  is  incident  on  the  end 
wall  and  gives  rise  to  a  reflected  wave.  The  potential  for  the  acoustic 
mode  is 


H  =  A  (1)  J 
1  P 
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o  p 


(1) 
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ik 


+  B 
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J  (7  (1)r) 

o\  p  I 


-ik  (1)(z-d) 
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(A.  56) 
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Near  the  end  wall  (z  =  d)  the  plane  viscous  and  thermal  waves  are  from 
(A.  12) 


(A.  57) 


The  velocity  and  temperature  can  be  calculated  just  outside  the  accom¬ 
modation  layer  using  (A.  4)  through  (A.  6)  to  get 
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(A.  58) 


The  boundary  conditions  on  r  and  T  are  the  same  as  those  given 
in  (A.  51).  However,  the  z  boundary  condition  must  be  modified  for 
acoustic  transmission  which  is  similar  to  (A.  18) 

z  =  ia_  wk,  7  ■  — ^  C  .  (A.  59) 

E  1  p  c 
w  w 

Here  a ^  is  a  constant  of  proportionality  which  depends  on  how  the  end 
wall  is  clamped  and  its  geometry.  One  of  the  problems  encountered  ex¬ 
perimentally  is  that  plate  resonances  in  the  end  walls  must  be  avoided. 
Otherwise  a  ^  is  not  a  constant  and  it  is  not  possible  to  separate  cavity 
resonances  from  end  wall  resonances.  As  before,  the  velocity  and 
temperature  in  the  gas  can  be  set  equal  to  the  velocity  and  temperature 
at  the  wall  to  give  a  set  of  equations  for  the  reflection  coefficients  r 

B  (i) 

p 

— =  r..  Thus  (A.  58)  becomes 

A  (1)  1 
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h  _a  uk  j  pc  \  k  (1)r  +k  (2)r  +k  (3)y  (1)2r 
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Solving  for  the  reflection  coefficient  by  Cramer's  rule,  we  get 


e  k  (1)-  h  ik  (1)d  ik,(1)d 

r  =  — ! -  e  P  =  e  P 

1  (1) 

e  k  1  '+h 
-  P 


where  and  €  are  given  by 


■±  =  1+  iaEkj 


P  C 
w  w 


(A.  60) 


(A.  61 ) 


(A.  62) 


Thus,  a  wave  of  amplitude  A^  originating  at  z  =  0  incident  on  an  end 
wall  located  at  z  =  d  gives  rise  to  a  reflected  wave  of  amplitude  r,  A^*) 
at  the  end  wall.  Notice  the  effect  of  traveling  a  distance  d  along  the  wave 
guide  has  been  in  the  reflection  coefficient.  In  addition,  k  (1),  e+,  and 
€  _  have  been  defined  for  algebraic  convenience. 
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The  field  in  the  cavity  is  sensed  at  the  opposite  wall  from  the 
source  wall  (z  =  d).  The  requirement  for  isolation  of  the  transducer 
from  the  cavity  leads  to  a  large  coupling  loss  between  the  input  signal 
and  the  output  signal.  Electrical  feedthrough  is  minimized  by  sensing 
the  field  in  the  cavity  at  z  =  d.  The  total  wave  traveling  toward  the 
end  wall  is 


ik  (1,d 
A  1  e  P 


ik  (1)d  i 2  k'  (1)d 


+  A  (1)e  P  e  P  +__  =  -£ 


m  ik  (1)d 

A  1  e  P 


ik'  {1)c 


.  (A.  63) 


1  -  2e 


Similarly,  the  sum  of  the  reflected  waves  at  z  =  d  is 


ik'  W, 


in  /,  v  ik'  (1,d 

(1)  ^P  d  A  (1)  i3kp  d  AP  6  P 

+  A  e  P  +  -  =  -2 - —  .  (A.  64) 
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The  total  excitation  at  z  =  d  is  the  sum  of  the  incident  and  reflected 
waves 
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(A.  65) 
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The  sum  of  two  angle  formulas  and  the  definition  of  k 
to  give  (note  e  +  e  =  1 ) 


,  (1) 
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can  be  used 
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The  final  expression  for  line  shape  of  a  cavity  is 


-G  J  <7  r) 
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(A.  67) 


where  we  have  used  (A.  56)  and  (A.  65)  to  simplify  (A.  66). 


The  working  relations  for  the  transport  property  determinations 
require  further  approximations  to  reduce  (A.  67)  to  a  manageable  form. 
The  general  magnitude  of  the  various  terms  indicates  how  such  approxi¬ 
mations  can  be  made.  When  the  real  part  of  kp(l)  d  =  niT  n  =  0,  1,  --- 
the  real  part  of  the  denominator  of  (A.  67)  becomes  very  small  and  the 
resonate  buildup  of  amplitude  is  observed.  The  radial  resonances  of  the 
cavity  occur  when  n  =  0.  That  is,  every  time  the  excitation  frequency  is 


( 


u 

,°P. 

c 


(A.  68) 


as  indicated  by  Eq.  (A.  45).  The  longitudinal  modes  occur  for  n  #  0  y  =0, 
so  that  P 


/  w 
|  nt> 

\  C 
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The  compound  modes  are  obtained  when  n  ^  0  and  y  $  0,  that  is,  the 
index  p  #  0  ^ 


(A.  69) 


The  approximate  signs  are  used  because  we  have  neglected  small  changes 
in  resonant  frequency  due  to  losses  in  the  cavity.  In  fact,  if  we  collected 
all  the  real  terms  of  the  denominator  in  Eq.  (A.  67)  and  set  them  equal  to 
zero  we  would  find  that  the  effect  of  the  first  order  correction  is  to  lower 
the  resonant  frequency  by  one -half  the  cavity  linewidth. 


The  wave  number  k  for  longitudinal  and  compound  modes  has  a 

large  k^^|  ~  1  cm  * j  real  part  and  a  small  ^1  m  <  10 

imaginary  part.  Therefore,  we  can  factor  out  k  ^from  the  denominator 
of  (A.  67)  to  give  ^ 


large 


(A.  70) 

The  transcendental  terms  in  the  denominator  of  (A.  70)  may  be  written 


sin  (kp^d)  =  sin  -  ~pjd  (A.  71") 

and 

coskp(1)d  =  (-l)n  cos(kp(1)-  “*  j  d  .  (A.  72) 

Near  longitudinal  or  radial  resonance  these  arguments  become  small  and 
Taylor  expansion  approximations  can  be  used. 

The  arguments  of  (A.  71)  and  (A.  72)  can  be  written  with  the  help 
of  (A.  46)  and  (A.  47) 
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near  the  resonant  frequency.  The  real  part  of  (A.  73)  is  of  the  order  of 


2  co 


6  -4  -3 

~10  ,  the  imaginary  part  is  of  the  order  of  10  .  Thus,  (A.  71) 


n 


and  (A.  72)  become 
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=  k  (1)d  +  0  (10'6) 
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(A.  74) 
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=  1  0  (10“6)  . 


(A.  75) 


The  line  shape  within  one  linewidth  of  the  center  frequency  reduces  to 
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(A.  76) 


The  linewidth  can  be  identified  by  using  (A.  73)  to  put  (A.  76)  into 
the  form 
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where  to  the  first  order 
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(A.  78) 


This  linewidth  is  the  frequency  difference  between  45  phase  shift  points. 
The  second  order  terms  are  again  lumped  into  multipliers  correcting  the 
viscosity  and  thermal  conductivity.  The  total  correction  terms  including 
those  originally  noted  in  Table  A.  6  are  given  in  Table  A.  8. 


*6  is  the  linewidth. 
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The  radial  mode  linewidth  required  a  slightly  different  approach 

because  k  ^ — *■ 0  fast  resonance  made  the  Taylor  expansion  (A.  73) 
p  2  2 

(1)  “  _2 

impossible.  The  wave  number  'd  goes  to  ~  ^ —  <10  for 


co 


01 


the  first  radial  mode.  The  small  angle  approximations  can  be  made  for 

sin  k  ^  d  and  cos  k  ^  d.  1 
P  P 

the  present  approximation  is 


sin  k  ^d  and  cos  k  ^  d.  The  line  shape  near  the  line  center  (A.  67)  in 
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The  linewidth  is  found  as  before  to  be 


.  1/2  .  to  p  c 

o  —  co  0  4-  a  7 d  — ^ — 

po  po  s  c  P.„c 


(A.  80) 
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The  second  order  terms  considered  as  corrections  to  viscosity  are  given 
in  Table  A.  8  for  the  radial  mode. 
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APPENDIX  B 


EFFECTS  OF  MULTIPLE  WAVES  IN  COUPLING  HOLES 


The  sound  is  coupled  into  the  cavity  through  small  holes.  The 
coupling  holes  are  themselves  wide  band  resonators.  That  is,  short 
cylindrical  wave  guides  terminated  on  one  side  by  the  electrostatic 
transducer  and  on  the  other  side  by  a  large  change  in  radius.  Both 
terminating  impedances  reflect  most  of  the  energy  incident  on  them. 
The  lowest  resonant  frequency  of  the  coupling  holes  is  (from  three  to 
six  times)  well  above  the  cavity  resonances.  In  addition,  because  of 
the  small  diameter  the  losses  are  50  times  greater  than  the  cavity 
losses.  The  Q  of  the  coupling  holes  is  at  most  0.  1  times  the  cavity,  Q. 


The  cavity  plus  transducers  consist  of  finite  cylindrical  wave 
guides  connected  together  as  shown  in  Fig.  Bl.  The  problem  is  greatly 
simplified  by  the  fact  that  all  the  transmission  coefficients  (T)  are  very 
small.  The  various  components  of  signals  reaching  the  receive  trans¬ 
ducer  can  be  classified  by  the  power  of  the  transmission  coefficients 
which  multiply  them.  Since  the  transmission  coefficients  are  always 
less  than  10“^  such  a  series  will  converge  fast. 


The  total  forward- going  wave  due  to  sound  initially  trapped  in 
the  transmitter  can  be  expressed  as  the  product  of  the  amplitude  of  the 
signal  introduced  at  z  =  0  (Aq)  times  the  response  of  the  coupling  hole  Fj 


A  =  A  F  = 
f  o  1 


A  e 
o 


-iw  (1>i, 

pi  1 


1 +RooRu  2i(kp(U<  i  +  *n) 


(B.  1) 


This  is  the  total  primary  wave  incident  on  the  interface  (1)  of  Fig.  Bl. 
The  expression  for  F^  is  slimming  multiple  reflections  of  the  waves 
trapped  in  the  coupling  hole  and  evaluating  at  z  =  l  .  The  reflection 
coefficients  Rqq  and  Rjj  are  defined  in  Fig.  Bl.  <p  ^  is  the  phase  shift 
of  the  reflection  coefficient  at  the  exit  of  the  coupling  hole  interface  (1). 
Similar  quantities  are  defined  for  each  of  the  three  parts  of  the  cavity 
resonator  in  Fig.  Bl. 

The  total  primary  wave  exiting  the  cavity  is  Af  Tj.  This  wave 
undergoes  multiple  reflections.  The  resulting  primary  forward  wave 
at  interface  (2)  is  then 


Bf  =  AoFlTl.F  • 


(B.  2) 
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Fig.  Bl.  Definition  of  terms  for  coupling  hole  problem. 
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This  signal  is  then  introduced  into  the  receive  coupling  hole  to  give  a 
total  received  primary  wave 


C,  =  A  F  T.  FT  ,F,  . 

*  O  1  It  i  b  2 


(B.  3) 


A  backward  primary  wave  is  also  generated  in  the  cavity.  This 
wave  incident  on  interface  (1)  is  transmitted  into  the  transmitter  coupling 
hole.  A  secondary  forward  wave  is  thus  generated  which  leads  to  a  sec¬ 
ondary  component  in  the  received  signal.  This  signal  has  two  additional 
transmission  coefficients  as  multipliers  since  it  passes  through  inter¬ 
face  (1)  two  additional  times  as  compared  with  the  primary  waves.  An 
additional  secondary  wave  is  generated  by  the  backward  primary  wave 
from  the  receive  coupling  hole.  The  total  secondary  wave  is 


R.  T 


1  2. 


F2F  exp 


-l(‘ 


(1) 


d  -  k 


(1) 


p2 


+  R0TlTl.rFl  ex»-i(kpl<1’,l-kp<1)d)JCf  (B'4) 

The  products  T2  T  ^  and  T  1  Tj  are  typically  of  the  order  of  10  ^  so 
that  secondary  and  higher  order  waves  can  be  neglected. 

The  receive  signal  is  given  by  Eq.  (B.  3).  The  reflection  and 
transmission  coefficients  are  complex  and  slightly  frequency  dependent. 
However,  this  effect  is  small  and  can  be  considered  as  a'  constant  sys¬ 
tem  phase  shift  over  a  cavity  linewidth.  The  reflection  and  transmission 
coefficients  at  the  cavity  coupling  hole  interface  are  the  traveling  wave 
coefficients  calculated  in  Appendix  C.  Notice  that  the  large  insertion 
loss  in  the  cavity,  typically  120  dB,  is  required  to  make  the  contributions 
of  secondary  and  higher  order  waves  small.  The  results  rapidly  become 
intractable  as  the  higher  order  terms  become  important. 
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APPENDIX  C 

REFLECTION  AND  TRANSMISSION  COEFFICIENTS  AT 

COUPLING  HOLES 


The  excitation  of  a  fluid  filled  cavity  with  coupling  holes  is  a 
complicated  process  because  small  holes,  the  order  of  the  wavelength 
of  viscous  and  thermal  waves,  must  be  used  to  preserve  the  Q  of  the 
cavity.  Each  of  the  three  wave  systems:  acoustic,  viscous  (vortex) 
and  thermal  (entropy)  waves,  satisfy  a  wave  equation  independently  at 
every  point  within  the  fluid.  The  waves  couple  on  interfaces  between 
solids  and  fluids.  This  coupling  dictates  the  relative  partition  of  energy 
between  the  acoustic,  viscous  and  thermal  wave  types.  However,  the 
acoustic  mode  wave  functions  alone  may  be  used  to  match  velocity  and 
pressure  at  the  interface  of  the  cavity  and  the  coupling  hole.  The  main 
contributions  of  viscous  and  thermal  waves  are  contained  in  the  imag¬ 
inary  parts  of  k  and  y  (1)  from  Appendix  A. 


In  the  present  discussion  we  calculate  the  linewidth  of  a  cavity 
fed  by  a  centered  coupling  tube  on  one  end  wall  and  tapped  by  a  centered 
receiving  tube  on  the  other  end  wall.  The  coupling  holes  cause  line 
broadening,  shift  the  center  frequency  and  introduce  a  phase  shift  be¬ 
tween  the  electrical  signals  and  the  pressure  wave  in  the  cavity.  The 
reflection  coefficient  for  traveling  waves  in  the  coupling  tube  and  the 
cavity  will  be  calculated  using  normalized  waves  from  Appendix  A, 
which  already  satisfy  the  boundary  conditions  on  the  side  walls  of  the 
cavity  and  the  coupling  tubes.  This  means  symmetric  radial  modes  for 
a  circular  wave  guide  in  the  present  case.  The  terminology  is  indicated 
in  Fig.  Cl. 


Here  we  consider  a  wave  incident  on  the  coupling  hole  from  the 
cavity  in  the  sth  radial  mode.  In  the  reflection  process  other  modes, 

( pn ,  are  generated  in  the  cavity  and  radial  mode  waves,  are  excited 

in  the  coupling  hole.  The  other  modes  in  addition  to  the  incident  mode 
are  required  because  the  particle  velocity  in  the  junction  plane  (CT )  be¬ 
tween  the  cavity  and  the  coupling  tube  does  not  correspond  to  a  pure  mode 
of  the  cavity  or  the  coupling  tube. 

We  insist  that  the  medium  does  not  come  apart  so  the  particle 
velocity  is  continuous,  and  that  energy  is  conserved  so  the  stress  is 
continuous  across  any  surface.  Losses  due  to  the  thermal  and  viscous 
waves  required  on  the  solid  portion  of  the  wall  (S)  are  taken  into  account 
by  adding  the  appropriate  linewidth  for  a  closed  cavity  to  the  linewidth 
due  to  losses  to  the  coupling  hole  area  (a)  (see  Fig.  Cl).  Note  that  in 
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Fig.  Cl.  Terminology  for  hole  coupling  problem. 
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what  follows  little  explicit  mention  of  thermal  and  viscous  waves  is 
made.  Since  no  coupling  of  the  viscous  thermal  and  acoustic  waves 
occurs  in  (7,  the  stresses  and  particle  velocity  due  to  each  wave  type 
separately  must  be  continuous  in  0 .  Viscous  and  thermal  waves  suf¬ 
ficient  to  satisfy  the  boundary  condition  on  S  and  the  side  walls  of  the 
tube  come  close  (1  pt  in  10^)  to  satisfying  the  continuity  condition  for 
all  the  waves  automatically.  The  effects  of  viscous  and  thermal  waves 
therefore  are  assumed  to  be  accounted  for  by  using  acoustic  traveling 
wave  modes  with  wall  losses.  Our  first  task  is  to  set  up  and  solve  an 
integral  equation  for  the  particle  velocity  uQs  in  (a)  when  the  sth  cavity 
mode  is  incident  on  the  cavity  end  wall.* 

Since  the  radial  part  of  the  traveling  wave  functions  for  tubes 
can  be  normalized  to  form  a  complete  orthogonal  set,  continuity  of 
axial  particle  velocity  (u0s)  at  z  =  0  (on  a)  can  be  assured  by  simply 
expanding  uQs  in  the  two  sets  of  traveling  wave  functions.  For  in¬ 
stance,  the  axial  particle  velocity  in  the  cavity  us  is 


-  j  K  z 
J  n 

uS  =  d>  sinh  (0  +  j  K  z)  +  V''  U  8  6  e 

Y  s  s  J  s  Lu  n  Y  n 

n*  s 


(C.  1) 


where  the  (sinh)  notation  is  used  to  represent  the  sum  of  the  incident 
and  reflected  sth  mode  wave.  Similarly,  the  axial  particle  velocity  in 
the  coupling  tube  is 


m 

u 


-2>, 

m 


0 

m 


jk__z 

m 


(C.  2) 


where  and  are  the  wave  numbers**  for  the  n  and  mth  modes  in 

the  cavity  and  the  coupling  tube,  respectively.  Note  that  we  have  assumed 

that  the  wave  equation  is  separable  with  the  solution  to  the  radial  wave 

equation  being  This  is  not  possible  in  the  comers  of  the  wave  guide. 

However,  a  small  error  discussed  below  is  introduced  (1  pt  in  10^)  by  this 

assumption.  The  coupling  coefficients  Uns  and  Ums  are  defined  so  that 

they  give  the  proper  particle  velocity  distribution  uQs  in  the  coupling  area 

<7.  That  is,  the  U  8  and  V„8  are  the  projections t  of  u  8  on  y  and  d) 

n  in  o  Til  1  n 


*The  notation  and  approach  are  those  of  Ref.  Cl. 

**The  wave  numbers  are  given  by  Eq.  (A.  46)  of  Appendix  A. 
notation  has  been  adapted  here  for  simplicity, 
t  The  abbreviated  notation  means 

<  ^  |  u  >  =  /  (^ud<7;dv  =  2TTrdv. 

JQ 


A  new 
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U 


and 


n 


m 


<  0  |  U  S  > 

r  n  o 


<  *  |  u  S> 

m  o 


(C.  3) 


(C.  4) 


The  generalized  phase  angle  (9  )  in  terms  of  the  coupling  coefficients  is 


sinh  9  =  <  0  u  >  . 

n  r  n  1  o 


(C.  5) 


Finally,  expressions  for  the  reflection  coefficient  Rg  for  the 
incident  wave  in  the  cavity  as  well  as  the  transmission  and  excitation 
coefficient  for  the  various  waves  excited  during  reflection  will  make 
the  above  definitions  more  concrete.  The  reflection  coefficient  is 


-29 


=  e 


(C.  6) 


The  excitation  coefficient  of  coupled  modes  in  the  cavity  is 


-29 

R  8  =  2U  8  e  n 
n  n 


(C.  7) 


and  the  transmission  coefficient  into  the  coupling  tube  is 


-0 

_  s  __  s  m 
T  =  V  e 
m  m 


(C.8) 


In  light  of  the  definition  of  9  n,  Un  and  Vm  above  all  the  needed  coef¬ 
ficients  can  be  calculated  once  the  particle  velocity  in  the  area  <T  is  known. 

Continuity  of  stress  reduces  to  continuity  of  pressure  for  the 
acoustic  wave.  The  normalization  implicit  in  Eqs.  (1)  and  (2)  means  that 
if  axial  velocity  due  to  the  mth  mode  is 


u  =  0  e 
m 


j  Kz 


(C.  9) 


then  the  pressure  due  to  the  mth  mode  is 


P  = 


0) 


v  0  e 

K  v  m 


j  Kz 


(C.  10) 
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The  pressure  in  the  cavity  is 


v  ^  s  U  -  j  K  z 

t  =  "  —  cosh»s+jKsz)-p2  n 

s  s  n 


(C.ll) 


Similarly,  in  the  coupling  tubes  the  pressure  is 


2 

w 


m 


J  k  z 

m  m 

^  e 
K  m 
m 


(C.  12) 


The  equality  of  the  pressure  at  z  ,  0  gives 

0s  cosh0  s  -  <G“  [uoS>  -  <G+(uoS>  =  0  (C.  1 3) 

where 

n^s  n 

G+=Eir"  *m(ro>  *n{T)- 
m  m 


(C.  14) 


(C.  15) 


and  rQ  is  an  independent  radial  coordinate  for  integration  purposes.  In 
order  to  put  (C.  13)  in  a  more  convenient  form,  we  multiply  both  sides  of 
(C.  13)  by  uQs  and  integrate  over  the  remaining  (r)  dependence  on  a.  In 
addition,  we  add  to  (C.  13)  after  integration  the  relation 

cosh©  <  u  S 1 0  >  =  sinh  20  (C.  16) 

SOS  u  s 

(obtained  from  (C.  5)).  The  final  integral  equation  for  u  8  is  then 


—  sinh  20 


=  2  cosh  0 


<  u  I  6  > 

o  r  s 


u 


>  +  <  u 


lG+|uoS>], 


(C.  17) 
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A  simple  variational  approach  to  (C.  18)  is  possible  because  the 
integral  equation  (C.  17)  is  quadratic  in  uqs.  A  useful  set  of  functions 
to  expand  uqs  for  the  variation  is 

30 

u  8  =  V'  A  v  cosh  0  (C.18) 

o  jji  n  s 

H=0 


where 


(C.  19) 


The  first  two  such  functions  are  used  in  the  present  work.  The  equations 
for  A^  obtained  from  (C.  18)  and  (C.  17)  as  described  below  are  symmetric 
and  easily  extended  to  any  order  of  p .  However,  calculation  for  p  =  0-2 
shows  no  significant  improvement  over  the  result  for  p=  0  and  1.  This 
assumed  functional  form  (C.  18)  can  be  put  into  (C.  17)  to  give 


1 

2 


tan  2d 

s 


2  A  <  1 1  >  +  2  A  <  X  ,  U  > 

o  1  1  s  1  *  1  r  s 


-Aq2  <  1 1  G’  +  G+ )  1  >  -  Aj2  <Xj|  G"  +  G+|Xl>  (C.  20) 


-  2  A  A  <  l|  G"  +  G+|x  >  • 

o  1  1 


As  explained  in  Ref.  1,  page  155,  integral  equations  such  as 
(C.  20)  can  be  used  to  obtain  the  "best"  values  of  the  unknown  parameters 
A  by  requiring  that  tan  2  6  s  have  a  minima  with  respect  to  the  variation 
or  the  A  .  Thus,  setting  the  partial  derivatives  with  respect  to  A0  and 
Aj  equal  to  zero  two  simultaneous  equations  for  A0  and  A^  are  obtained 
(note  we  have  used  the  symmetry  of  G+  and  G  ) 

<l|G_+G+|l>  Aq+<1  |G_+G+|x1>A1  =  <'lUfl> 

(C.21) 

<l|  G"  +  G+|Xl  >  Aq+  <  xxl  G"  +  G+  |  Xl>  Ax  =  <X1kB>  • 


*See  Ref.  Cl  for  formal  results. 
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These  equations  can  be  solved  for  A 


1 


and  A 

o 


to  give 


Ao  =  f  [<1|^S>  <x1Ig-  +  g+Ix1>-<x1Usxi|g'  +  g+Ix1>] 

A1  =  1  I  G"  +  G+l  I  G“  +  G+|  Xl>]  (C.  22) 

where 

A  =  <  1  |  G"  +  G+ 1  1  >  <Xj  I  G"  +  G+l  Xt>  -  <l|  G'  +G+I  Xx>2  • 

Once  the  wave  functions  0  n  and  0  m  are  known  the  constants  AQ  and  Aj 
can  be  evaluated  from  (C.  22)  and  2  0  calculated  from  (C.  20).  This 
gives  us  the  required  reflection  coefficient  from  (C.  6).  In  addition, 
substituting  the  now  known  particle  velocity  u0s  from  (C.  2)  into  (C.  3) 
and  (C.  4)  all  the  coupling  coefficients  Uns  and  Vns,  and  the  excitation 
(C.  7)  and  transmission  coefficients  can  be  calculated. 


The  wave  functions  required  to  evaluate  the  integrals  indicated 

(1). 


above  are  the  traveling  wave  functions  Jq  (y^ 


f)  of  Appendix  A.  These 

wave  functions  must  be  made  orthogonal^)  and  normalized  before  they 
can  be  used.  Orthogonalization  of  the  cavity  wave  functions  is  not 
important  and 


*P  ■  T 

r  a 


JoV> 


However,  the  suitably  normalized  wave  functions  for  the  coupling  tube 
functions  given  in  Table  C.  1  are  considerably  more  complicated. 


Fortunately,  only  the  first  two  functions  0  are  necessary.  Re¬ 
calling  the  size  of  y  from  Appendix  A  (Eq.  A.  47)  the  orthogonalized 
wave  functions  are  seen  to  reduce  to  the  lossless  wave  functions  as 
2 

y  ^  — ►  0.  In  fact,  by  the  time  p  =  100  atm  the  lossless  wave  functions 
describe  the  reflection  and  transmission  coefficients  typically  to  1  part 
in  105. 
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Table  C.  1.  First  two  orthogonalized  wave  functions  in  a 
coupling  tube. 
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